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1. 


INTRODUCTION 


It  is  known  that  the  surface  currents  associated  with  travelling  in¬ 
ternal  waves  can  be  sufficiently  intense  to  produce  modifications  of  the 
surface  wave  structure  that  are  visible  both  to  the  eye  and  to  radar 
(Ziegenbein,  1970;  Perry  and  Schimke,  1965;  Polvani,  1972).  In  a  series 
of  measurements  off  San  Diego,  Lafond  (1963;  1966)  has  shown  that  visible 
surface  slicks,  under  certain  wind  conditions,  move  with  the  underlying 
train  of  internal  waves.  A  later  series  of  measurements  at  the  same  in¬ 
strument  platform  has  shown  that  these  slicks  are  detectable  with  radar 
and  exhibit  phase  correlation  with  the  current  pattern  produced  by  the 
internal  wave  train  (Polvani,  1972). 

Two  different  mechanisms  have  been  suggested  as  responsible  for  the 
coupling  of  internal  waves  to  surface  waves,  one  associated  with  changes 
in  surface  composition  and  the  other  a  direct  modification  of  wave  struc¬ 
ture.  Lafond  (1966)  suggested  that  the  divergence  and  convergence  of  the 
surface  current  over  the  leading  and  trailing  portion  of  an  internal  wave 
maximum  disperses  or  concentrates  extant  organic  material  on  tne  surface 
in  a  way  correlated  with  the  internal  wave  phase.  The  resulting  increase 
in  the  damping  coefficient  in  the  convergence  zones  reduces  the  amplitude 
of  short  wavelength  gravity-capillary  waves  in  these  areas  and  changes  the 
average  wave  slope  and  reflectivity. 

A  number  of  authors  have  noted  the  resonant  effect  that  may  be  ex¬ 
pected  on  a  surface  wave  packet  that  travels  at  the  same  speed  as  the  in¬ 
ternal  wave  train  (Hartle  and  Zachariason,  1969;  Holliday,  1971;  Rosen- 
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bluth,  197J ;  Phillips,  1971;  Ko,  1972).  Since,  at  resonance,  the  wave  packet 
is  able  to  tra  -l  great  distances  while  remaining  at  the  same  phase  point  of 
the  internal  wave,  the  continuing  interaction  leads  to  a  continuous  energy 
transfer  from  the  internal  to  the  surface  wave. 


The  theoretical  models  of  the  wave-current  interaction  may  be  divided 
into  two  broad  categories,  which  we  might  characterize  as  space  oriented 
(wave  packets)  or  wavenumber  oriented  (modal). 


The  wave  packet  modpls  (WKB,  eikonal,  ray  theory)  developed  by  Longuet- 
Higgins/Stewart  (1961),  Whitham  (1966),  Phillips  (1966),  Holliday  (1971), 

Ko  (1972),  consider  the  surface  as  a  superposition  of  a  number  of  spatially- 


localized  wave  packets,  each  characterized  by  a  characteristic  wavenumber 
k(x,t)  and  wave  frequency  u(x,t).  These  packets  move  along  trajectories 
defined  by 


dx 

dt 

dk 

dt 

where  the  notation  indicates 


3to 

■+ 

9k 

.  _  3w 

9x 

gradients  in  k  and 


x-space 


respectively. 


The  mode-oriented  models  [Ball  (1964),  Hasselman  (1967),  Hartle/ 
Zachariason  (1969),  Rosenbluth  (1971))  describe  the  surface  essentially  as 
a  superposition  of  various  modes  (sine  waves  for  the  unperturbed  ocean) 
and  concentrate  their  attention  on  the  transfer  of  energy  between  modes. 
When  applied  to  small  amplitude  waves,  both  treatments  are  equivalent  since 
the  requirements  for  validity  of  the  WKB  approximation  are  usually  well 
satisfied  for  the  waves  of  interest.  However,  the  equivalence  tends 
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to  be  obscured  in  development  of  computational  models.  In  a  modal  descrip¬ 
tion,  spatial  structure  is  bound  up  in  the  details  of  the  phase  relation¬ 
ships  between  different  wavenumber  modes.  On  the  other  hand,  the  wave¬ 
number  structure  and  spectrum  is  related  to  the  details  of  the  spatial 
correlation  function  in  the  wave  packet  description.  Rosenbluth  (1971) 
combines  elements  of  both  descriptions  in  a  "full-wave"  treatment  of  a 
linear  wave  interacting  with  a  variable  surface  current  and  shows  the 
equivalence  between  eigenmodes  and  wavepackets  trapped  over  troughs  of 
the  internal  wave. 

A  major  difficulty  arises  in  the  application  of  these  mechanisms  to 
describe  ocean  wave  spectra.  The  ocean  surface  under  typical  wind  condi¬ 
tions  is  poorly  described  as  a  simple  linear  superposition  of  sine  waves. 

On  the  contrary,  the  wave  amplitudes  and  slopes  are  usually  driven  to  the 
point  where  non-linearities  are  dominant  in  providing  a  saturation  limiting 
of  the  wave  heights.  These  effects  manifest  themselves  in  the  appearance 
of  sharp  crested  waves  and  wave  breaking  (white  caps).  In  the  present 
series  of  reports  we  will  develop  a  description  of  the  dynamical  properties 
of  ocean  wave  spectra  involving  a  transport  or  relaxation  equation.  Here 
we  write  the  rate  of  change  of  the  mean  or  expected  value  of  the  wave  height 
power  spectrum  in  the  form 

—  F(x,k,t)  =  A(x,k)  V+U  +  generation  terms  -  loss  terms. 

at  x 

where  F  is  assumed  to  depend  on  (x,k,t)  and  U  is  the  local  time  averaged 

x 

surface  current  gradient.  The  detailed  development  and  modeling  of  the  var- 
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ious  phenomena  contributing  to  wave  growth  and  decay  will  be  discussed  in 
subsequent  reports.  In  the  present  report  we  wish  to  lay  the  groundwork 
for  the  discussion  of  the  wave-current  interaction  and  to  review  the  physics 
of  the  interaction  process.  To  do  this  we  will  discuss  primarily  a  linear 
description  of  the  ocean  surface.  Many  of  the  non-linear  effects  coupling 
various  waves  together  can  be  considered  as  interactions  to  various  orders 
of  different  wave  fields  and  can  be  understood  within  the  context  of  a 
linear  model  in  terms  of  the  interaction  of  a  linear  wave  with  a  prescribed 
current  field.  Thus,  part  of  the  motivation  of  the  present  work  is  to 
describe  the  resonant  interaction  of  a  surface  wave  with  a  moving  current 

pattern  and  part  to  clarify  some  of  the  concepts  that  underlie  the  non¬ 
linear  theory. 

We  consider  the  current  produced  near  the  ocean  surface  by  a  long 
wavelength  internal  wave  propagating  along  a  sharp  thermocline  at  depth  D 
(see  Figure  1).  Fcr  small  percentage  changes  in  the  density  across  the 
interface  (A p/p),  the  currents  and  potential  of  the  internal  wave  are 


»(0) 

A 

_  0 

K 

cosh  Ky 

sin  K(x  -  cQt)  , 

(1.1) 

j(°) 

X 

=  A 

0 

cosh  Ky 

cos  K(x  -cQt)  , 

(1.2a) 

,(0) 

y 

-  A 

0 

sinh  Ky 

sin  K(x  -  c  t)  , 

(1.2b) 

where  cq  is  the  phase  velocity  of  the  interfacial  wave.  The  vertical 
velocity  of  the  surface  has  been  taken  to  be  zero  (valid  in  the  limit 
Ap/p  0).  In  this  paper  we  study  the  effect  of  this  current  [Equation 

(1.1)]  on  the  motion  and  pattern  of  the  small  wavelength  surface  gravity 
waves . 
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To  get  a  feel  for  the  basic  physical  versus  the  mathematical  aspects 
involved  in  the  problem,  we  examine  first  the  dispersion  relation  for  the 

surface  waves  for  a  uniform  (deep) -ocean  in  the  presence  of  a  uniform  cur¬ 

rent.  It  is  well  known  that  the  dispersion  relation  in  a  coordinate  sys¬ 
tem  in  which  the  ocean  is  stationary,  has  the  form 

“2  *  S*  .  (1.3) 

or 

to  =  ±i/gk 

where  u>  is  the  angular  wave  frequency,  k  the  wavenumber  of  the  surface  wave 
and  g  the  gravitational  acceleration.  The  velocity  of  a  point  of  constant 
phase  on  the  wave,  i.e.,  phase  velocity,  is 


c  =  &  =  ± /g/k  . 

ph  a) 


(1.4) 


This  dispersion  relation  is  plotted  in  Figure  2a.  Obviously,  for  each  value  of 
|u|  there  are  only  two  possible  solutions  and  two  wave  velocities  (±  g/|o»|)  cor¬ 
responding  to  waves  travelling  to  the  right  or  to  the  left. 

If  we  now  make  a  Galilean  transformation  to  an  observer  moving  over 
the  ocean  surface  at  velocity  cq,  events  appear  to  have  a  different  time 
dependence.  We  must  replace  oi  by  u>  +  cQk  where  cq  is  the  observer  velocity 
relative  to  the  ocean.  This  has  the  tendency  to  twist  the  dispersion  curves 
downward,  i.e.,  the  trajectory  of  the  ocean  surface,  which  was  the  real  axis 
u)<=0,  now  becomes  (see  Fig.  2b) 


-  c  k 
o 


(1.5) 


This  moving  observer  finds  himself  with  a  rather  more  complicated  description 
of  the  same  simple  wave  motion  seen  by  the  stationary  observer.  First  he 
finds  that  there  is  a  maximum  frequency  (u>*)  above  which  waves  will  not  propagate 
freely.  At  any  lower  frequency  he  finds  that  there  are  now  four  different  wave- 
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numbers  and  wave  velocities  that  will  propagate  freely  fcr  each  given  value  of 
|w|  Instead  of  two  as  in  the  ocean  fixed  frame.  If  this  observer  knew  that  in 
his  coordinate  system,  the  velocity  of  the  ocean  were  locally  uniform  and  inde¬ 
pendent  of  time,  then  he  would  correctly  conclude  that  possible  solutions  have 
the  form  exp[-iwt  +  ikx]  where  there  are  two  permissible  right  travelling  waves 
(w/k  >  0)  and  two  left-travelling  waves  (cu/k  <  0)  that  are  consistent  with  the 
same  value  of  |w|  . 

All  this  Ls  of  little  interest  if  the  ocean  current  is  truly  independent 
of  position  since  in  that  case  the  waves  are  simply  manifestations  of  the  two 
simple  surface  waves  in  Equ.  (1.3).  These  two  waves  (or  four  depending  on  the 
coordinate  system)  are  uncoupled  and  the  transformation  to  a  moving  coordinate 
system  merely  complicates  a  simple  problem.  However,  when  the  surface  current 
varies  with  position  (but  not  so  rapidly  as  to  invalidate  the  local  validity 
of  the  dispersion  relations)  the  problem  is  distinctly  different.  In  the 
particular  case  of  a  travelling  current  disturbance  there  is  a  unique  coor¬ 
dinate  system  in  which  the  surface  current  distribution  is  Independent  of 
time.  In  this  coordinate  system,  the  solutions  are  harmonic  in  time  (e  ^w^) 
but  the  spatial  gradients  may  permit  or  induce  coupling  of  the  four  waves. 

Of  particular  interest  are  the  two  right  travelling  waves.  Reference  to 
Fig.  2b  shows  that,  although  the  phase  velocities  for  both  are  in  the  direction 
of  the  internal  wave,  the  group  velocities  (du/dk)  have  opposite  signs. 

These  two  waves  lie  on  either  side  of  a  turning  point  (where  the  group  velocity 
■t  vanishes) .  It  will  be  shown  in  section  4  that  a  wave  packet  incident  on  a  finite 

region  of  steady  surface  currents  can  be  totally  reflected  for  certain  current 
patterns  and  that  the  two  right  travelling  waves  represent  the  incident  and 
reflected  wave  trains.  The  observer  in  the  ocean  coordinate  system  would 
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Interpret  this  reflection  as  follows.  A  surface  gravity  wave  packet  traveling 

at  velocity  c  (group)  overtakes  a  traveling  surface  current  region  [translat- 
8 

ing  at  velocity  c(  <c  )].  There  Vill  be  a  range  of  values  for  c  -c,  dependent 

6  ® 

on  the  maximum  amplitude  of  the  surface  current,  for  which  the  gravity  wave 
loses  energy  to  the  surface  current  to  the  extent  that  its  propagation  velo¬ 
city  is  reduced  to  2c  -  c  [i.e.,  a  velocity  decrease  of  2(c  -c)  because  of 
reflection].  In  the  same  manner,  waves  having  group  velocities  less  than 
but  sufficiently  close  to  surface  current  propagation  velocity  will  be  ac¬ 
celerated  by  the  front  edge  of  the  moving  current  field  to  a  velocity 

c +  (c  -  c  ).  Thus,  the  surface  current  field  acts  as  a  piston  to  a  certain 
8 

range  of  wavelengths  and  pushes  these  waves  along  ahead  of  the  current  field. 

Waves  near  the  condition  c  =  c  will  tend  to  pile  up  in  front  of  or  behind  the 

g 

reflection  point  and  form  caustics. 

This  type  of  interaction  may  be  thought  of  as  a  resonance  phenomenon 
where  waves  in  a  restricted  range  of  frequencies  interact  strongly  with  the 
external  forcing  function.  Indeed,  in  A^modal  description  where  the  surface 
is  represented  as  an  assembly  of  harmonic  oscillators,  the  reflection  phenomena 
described  above  manifests  itself  specifically  as  a  classical  resonance  inter¬ 
action,  between  the  oscillators  (surface  waves)  and  the  external  forcing  func¬ 
tion  (internal  wave). 

In  recent  years  considerable  interest  has  been  shown  in  this  type  of 
interaction  and  a  number  of  treatments  are  available  [Longuet-Higgins  and 
Stewart  (1960);  Longuet-Higgins  (1963);.  Ball  (1964);  Phillips  (1966); 
Hasselman  (1967);  Hartle  and  Zachariason  (1969);  Rosenbluth  (1972); 

Phillips  (1972)].  The  problem  also  has  in  interpretation  in  terms 
of  parametric  oscillator  theory.  In  spite  of  this  extensive  analysis, 
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the  relationships  between  the  various  methods  are  not  readily  transparent. 
Methods  based  on  the  WKB  approximation  allow  trajectories  of  individual 
packets  of  energy  to  be  followed.  In  available  treatments,  turning 
points  can  appear  in  the  solutions  where  the  WKB  approximation  falls  (blockage 
effects,  according  to  Phillips).  Modal  methods  based  on  the  Born  approxi¬ 
mation  show  a  related  singularity.  The  turning  point  here  appears  as  an 
oscillator  driven  at  resonance. 

It  is  the  purpose  of  this  paper  to  analyze  in  detail  the  specific 
case  of  linear  surface  waves  travelling  on  a  spatially  variable  but  steady 
current.  The  analysis  is  exact  and,  as  a  consequence,  permits  us  to  obtain 
bounded  solutions  in  the  neighborhood  of  the  resonance  regiun,  By  taking 
limiting  cases  of  these  exact  solutions,  we  are  able  to  demonstrate  the 
relationships  between  the  various  more  approximate  models  and,  in  particular, 
we  are  able  to  define  their  respective  ranges  of  validity. 

The  details  of  these  approaches  are  reviewed  in  the  following  section, 
along  with  an  introduction  to  the  present  approach  to  the  problem.  In  all 
case8»  detailed  theoretical  analysis  hast  been  limited  to  a  linear 
description  of  the  interactions.  Since  the  sea  surface  is  usually  in  a 
strongly  nonlinear  state,  such  approaches  may  have  limited  quantitative 
applicability.  However,  the  linear  analysis  provides  the  insight  necessary 
to  the  development  of  a  non-linear  theory. 

The  following  sections  segment  the  problem  of  internal  and  surface 
wave  interactions  as  follows:  Section  2  reviews  some  of  the  previous  work 
done  on  the  problem  to  familiarize  the  reader  with  the  difficulties  encountered 
and  to  point  out  where  these  studies  left  off.  Section  3  provides  an 
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introduction  and  review  of  scattering  theory  as  it  will  be  applied  to  the 
present  problem.  A  set  of  coupled  differential  equations  is  constructed 
from  an  application  of  this  theory  under  the  assumption  of  a  weakly  re¬ 
flecting  current.  These  equations  are  then  solved  analytically  for  a 
linear  current.  Section  4  discusses  solutions  in  both  the  Born  and  eikonal 
(WKB)  approximations.  These  approximate  solutions  are  then  related  f.o  the 
results  of  Section  J.  Section  5  summarizes  the  paper  stressing  the  inter¬ 
connection  between  the  differing  methods  of  approach  and  gives  a  physical 
interpretation  to  these  methods.  The  connection  with  previous  results  and 
the  interpretation  of  surface  slicks  in  terms  of  reflected  waves  is  also 
made. 
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2.  COUPLED  INTERNAL  AND  EXTERNAL  WAVES  (REV/ EW) 


We  consider  the  wavn  motion  at  the  surface  of  the  ocean  under  the  influ¬ 
ence  of  a  non-uniform  surface  current.  In  the  analysis  the  ocean  is  assumed  to 

he  homogeneous  and  irrotat lo,.,l .  Within  this  general  framework,  we  may  writ,  the 
momentum  equation  for  a  fluid  as 


du  1  _  -> 
dF  ~  •  p  7p  +  8 


(2.1) 


where,  under  the  assumption  of  irrotational  flow.  V  «  J  .  0.  The  juantities  in 

Eq.  (2.1)  are  defined  as  follows:  „  is  the  fluid  density.  £  is  the  gravitational 

acceleration,  p  is  the  pressure,  and  2  is  the  velocity  of  the  fluid.  If  write 

^  -  »♦.  where  ♦  is  the  potential  description  of  the  veiocity  field,  then  Eq.  (2.1 
may  be  written  as 


if  +  ;.v;  .  +  .V£tgyj 


(2.2) 


where  the  operator  d/dt  in  Eq.  (2.1)  has  been  replaced  by  the  Eulerian  derivative 
3/3t+u-V  Eq.  (2.2)  and  the  density  has  been  assumed  to  be  constant. 

Equation  (2.2)  may  be  integrated  immediately  to  yield  the  general  equation 


(2.3) 


*t  +  7  (7$)2  +  p  +  gy  =  f(t) 

where  in  the  following  we  set  f(t)  ■  n  *  n, 

B  e  sec  ru;  -  0.  At  the  surrace  of  the  fluid  we  will 

assume  that  the  pressure  is  zero  [p-0v-v  *  *  , 

ip  u,  y  m  ygt  $  -  *g],  so  that  Eq#  (2  3) 

becomes 


3$ 

Jt~  +  j  (7#)s  + 


0  . 


*  'it  COnditl°n  -  *****  that  hydrostatic 


(2.4) 


equi- 
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ThiH  equation  may  be  further  modified  by  recalling  that  the  rate  of  increase 
of  yg  following  a  fluid  element  is  the  vertical  component  of  the  fluid 
velocity: 


and  that 


dj> 

dt 


—  +  7$  .  7$ 


(2.6) 


Introducing  Equations  (2.5)  and  (2.6)  into  (2.4)  yields  the  differential 
equation 

{(ft  +  '  7)(fc  +  I  ”*  ■  v)*  +  8  !y  }s  ■  0  ■  <2'7> 

To  solve  Equation  (2.7)  we  expand  the  velocity  potential  as  follows: 

*  =  $(o)  +  *  ,  (2.8) 

where  4>  is  the  potential  associated  with  the  short  wavelength  surface  waves 
and 

Co )  A0 

v  -  —cosh  Ky  sin  K(x-ct)  .  (2.9) 

The  separation  given  by  Equation  (2.8)  essentially  couples  the  short  sur¬ 
face  waves  to  the  much  longer  internal  waves  by  superimposing  the  field  po¬ 
tentials. 


We  substitute  the  potential  from  Equation  (2.8)  into  Equation  (2.7) 
and,  retaining  only  up  to  the  linear  terms  in  obtain 
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[f^  +  It  (V$°*7)*  +  I7  (^*7)*°  +  74°. V  |i  +  (74°  *7)  (74>° *7)4> 


+  7$°  •  7 (7<p  *7)$°  +  7^*7  |~ 


+  ^•v(^.y,t-|aurface  +  (8  |i) 


=  0 


surface 


(2.10) 


This  expression  may  be  Immediately  reduced  by  using  the  condition  A  <<  c  and 

o 

noting  that  the  wavenumber  k  for  the  surface  wave  is  very  much  greater  than 
the  wavenumber  K  for  the  internal  waves,  i.e.,  K  <<  k.  With  these  two  con¬ 
ditions,  Equation  (2.10)  becomes 

[~2  +  8  If  +27*(0)*7|£l  -  0  (2.11) 

L  -» surface 

which,  together  with  the  Laplaclan 


0 


(2.12) 


defines  the  surface  wave  potential  for  our  problem.  Note  that  we  have 
assumed  the  wave  to  be  uniform  in  the  z-direction. 


Let  us  carry  out  the  remaining  analysis  in  terms  of  an  arbitrary 
traveling  internal  wave  packet  74>°  -  U(x-ct).  For  internal  waves  whose 
wavelength  is  long  compared  to  the  depth  of  the  interface  along  which  they 
propagate,  the  waves  are  non-dispersive  and  arbitrary  functions  U(£) , 

[£  =  x-ct],  may  be  considered.  This  function  is  analogous  to  the  potential 
produced  by  a  massive  particle  from  which  much  lighter  particles  (surface 
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waves)  are  being  scattered.  Reference  to  this  picture  will  make  the  con¬ 
cepts  introduced  in  the  following  analysis  more  transparent. 

Ball  (1964)  showed  that,  for  a  two-fluid  model  of  the  ocean,  resonant 
interactions  involving  two  surface  waves  and  one  internal  wave  can  occur. 
His  argument  was  basically  that  non-linear  quadratic  terms  in  the  wave  equa¬ 
tion  can  cause  a  combination  wave  (frequency  u>2 ,  wavenumber  k2)  to  be  gen¬ 
erated  by  the  beating  of  two  simple  harmonic  waves: 

o>2  =  ±  m 

and 

.  ,  x  (2.13) 

k2  »  k^  ±  K 

If  this  combination  [Equation  (2.13)]  wave  matches  one  of  the  natural  modes 
of  the  unperturbed  ocean,  then  a  resonance  condition  exists. 


To  follow  this  argument,  we  consider  the  wave  equation  valid  for  low 
speed,  low  amplitude  flow  [Equation  (2.11) J,  which  we  can  write*  as 


LA  +  R  1± 

o  '  K  * 


at 


ay 


-  2U  <*■'>  f^sF  • 


(2.14) 


When  the  right  hand  side  of  Equation  (2.14)  is  approximated  by  the  unperturbed 
surface  wave  solution  (the  Born  Approximation) 


$0  ~  a  cos  (kj^x  -u^t) 


(2.15) 


*  U  (x.t)  is  assumed  to  vary  slowly  over  distances  and  times  in  which  <f> 
oscillates.  Thus,  time  and  space  gradients  of  U  (x,t)  are  neglected. 
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and  the  travelling  current  field  by 


U(x,t)  =  A  cos(Kx  -  tit ) 


(2.16) 


[where  we  have  used  Equation  (2.9)],  combination  waves  are  generated.  We 
may  use  Equations  (2.15)  and  (2.16)  in  Equation  (2.14)  to  obtain 
2 

^^2  +  *  3y"  “  |cos|(k^+K)x  -  (oj^ +tj)tJ  +  cos£(k^-K)x  -  (o)^-o))t 

where  u  =  Kc.  (2*17) 


Of  particular  interest  is  the  long  wavelength,  slowly  varying  difference 
frequency  term.  Two  waves  having  almost  the  same  frequency  and  wavenumber 
and  travelling  in  the  same  direction  will  beat  and  give  rise  to  an  envelope 
that  has  a  much  lower  frequency  and  wavenumber.  Any  nonlinear  terms  will 
result  in  a  generation  of  a  wave  component  at  the  beat  frequency.  The 
wavenumber  of  the  beat  envelope  is  proportional  to  the  difference  between 
the  wavenumbers  of  the  two  primary  waves.  The  phase  velocity  of  the  beat 
wave  is  proportional  to  the  difference  between  the  frequencies  of  the 
primary  waves  and  can  be  shown  to  be  equal  to  the  harmonic  mean  of  the 
group  velocities  of  the  primaries. 

When  this  beat  wave  matches  the  internal  wave  in  frequency,  wavenumber, 
and  phase  velocity,  a  resonance  occurs  (see  below).  In  general,  the  coupling 
results  frcn  various  non-linear  terms.  Quadratic  terms  involving  two  sur¬ 
face  waves  can  create  combination  or  beat  waves  which  can  resonate  with  the 
internal  wave.  As  just  discussed  above,  quadratic  terms  involving  one  sur¬ 
face  wave  and  internal  wave  can  cause  coupling  to  a  second  surface  wave. 
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Marti e  and  Zacharlason  (1969)  gave  a  detailed  analysis  of  this 
problem;  that  is  the  solution  of  Equation  (2.14)  in  the  Born  Approximation, 
for  waves  travelling  in  the  same  direction.  The  right  hand  side  of  Equa¬ 
tion  (2.14)  was  evaluated  with  <j>Q  taken  as  the  unperturbed  surface  wave  po¬ 
tential  and  the  general  problem  of  an  arbitrary  travelling  current  wave  was 
treated.  They  also  gave  detailed  attention  to  the  variation  of  wave  height 
and  wave  slope  resulting  from  the  interactions. 

The  basic  features  of  their  solution  can  be  exemplified  by  treating 
the  special  case  of  a  sinusoidal  surface  current.  Equation  (2.17).  Since 
we  are  only  concerned  with  the  near  resonant  case,  we  consider  only  the 
difference  frequency  term,  replace  the  cosine  by  an  exponential,  and  find  a 
solution  to  the  equation 

+  g  =  aAk1u,1  exp^i^-lOx  -  (u^-Kc)^]  .  (2.18) 


We  may  write  <p(x,y,  t)  in  the  form 

<Hx,y,t)  =  ^  B(k)  exp  j ky  -  i  £kx  -  ojtJ  |  (2.19) 

k 

so  that  <p  +0  as  y+-  and  we  note  that  for  a  right  travelling  wave, 


M.  =  i  i! 

9y  "  9x  • 


The  solution  of  equ.  (2.18)  then  takes  the  form,  in  general, 


r  6 (u  -  u).  )6 (k  -  k.) 

♦  “  /  ,  J  a j  - > - L  du>  dk  , 

j 


u)  -  gk 


(2.20) 


(2.21) 
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where  u>j  and  k^  are  the  frequencies  and  wavenumbers  of  the  combination 

waves  generated  by  the  non-linear  terms,  i.e.,  k^-K,  and  uij— Kc  for  equ.(2.18), 

and  the  are  the  appropriate  coefficients.  Resonance  occurs  when,  for  some  j, 


•i  ■  8kj  ' 

In  our  example,  resonance  curs  if 


(2,2.2) 


(u>1  -  Kc;  =  g(kx  -  K) 


(2.23) 


Using  the  fact  that  =  gk^,  we  may  reduce  Eq.  (2.23)  to  the  form 

2 


2^-1  -  * 


(2.24) 


In  terms  of  the  surface  wave  group  velocity  c  (■  g/2u1),  we  may  write 

Si  x 


1  + 


-81 


4kl  (cgJ 


(2.25) 


or,  for  K/k^  <<  1, 


c 

C8l 


-  1 


(2.26) 


Thus,  resonance  obtains  when  the  group  velocity  of  the  surface  waves  is 
approximately  equal  to  the  internal  wave  velocity.  Note  that  this  state¬ 
ment  is  identical  to  the  statement  that  the  difference  wave  between  two 
similar  surface  waves  (which  could  be  generated  by  non-linear  effects) 
match  the  internal  wave  in  frequency,  wavenumber  (and  phase  velocity). 
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These  concepts  are  more  clearly  Illustrated  by  reference  to  the  (k,w) 

*► 

plane.  A  point  11  =  (k,aj)  in  this  plane  represents  a  general  wave.  The 
curve  which  is  the  locus  of  freely-propagating  waves  is  the  dispersion  re¬ 
lation  w  =  f (k) .  In  general,  quadratic  non-linear  terms  will  generate 
sum  and  difference  waves.  In  particular,  a  free  propagating  wave  defined 
by  the  vector  fi2[  =  (k^)]  can  couple  with  an  internal  wave  R[  =  (K,R)]  to 
produce  the  beat  wave 

$1  =  fi2  +  fl  . 

This  coupling  will  be  resonant  if  ^  is  a  characteristic  mode  of  the  surface. 
Three  such  combinations  are  shown  in  Figure  2c.  For  two  of  these 

•f 

the  combination  wave  (fi^)  does  not  fall  on  the  locus  of  freely-propagating 
waves.  However,  one  such  combination  does  satisfy  the  resonance  con¬ 
dition.  It  is  obvious  from  the  geometry  of  Figure  2c  that,  when  the  inter¬ 
nal  wave  has  low  frequency  (R  <<  R2>  K  «  k2>,  the  slope  R/K  (which  is  the 
internal  wave  phase  velocity  c)  is  essentially  equal  to  the  slope  of  the 
dispersion  curve  for  waves  satisfying  the  resonance  condition,  i.e., 

„  R  du) 

C  =  —  -  3  C 

K  dk  cg  *. 

To  examine  the  resonant  and  near  resonant  case  it  is  appropriate  to 
treat  an  initial  value  problem  where  we  assume  the  interaction  to  be 
turned  on"  at  a  certain  time.  For  the  interaction  starting  at  time  t=0, 

£  X  / 

i.e.,  <^(x,0)  *  ^  (x,0)  =  0,  the  solution  has  the  form,  [where  #  + 

a  expUkjX  -  iu^t)  ] , 
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Figure  2c,  surface  wave-internal  wave  interactions 
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AkjW^  a 


♦  ' 


i 


0)*  (to*  +  to) 


|  sin  to*t 


ito* 

to*  -  to 


1 


i(k1  -  K)x 

(2.27) 


where  u>  is  the  difference  between  the  surface  wave  frequency  and  the  in¬ 
ternal  wave  frequency  (to  *  to^  -  Kc)  and  to*  is  the  frequency  of  a  surface 
wave  having  a  wavenumber  equal  to  the  difference  between  that  of  the  sur¬ 
face  wave  and  of  the  internal  wave  Ito*  =  /g(k^  -  K)  J  . 


When  the  exciting  frequency  to  is  made  equal  to  the  frequency  to*,  ex¬ 
pand  Equation  (2.27)  and  we  find  that  the  wave  amplitude  resonates  and 
grows  linearly  in  time: 

k.to.  Aa  1 

(fr'Cto-to*)  =  i  1  ~--x —  7sin(io*t)  -  to*t  e 

2to*  I 


i(k.  -  K)x 


(2.28) 


Close  to  but  not  at  resonance,  the  root  mean  square  value  of  the  wave  po 
tential  is 


k^to^  Aa 
2oo*  (to*  +<o) 


(2.29) 


Hartle  and  Zachariason  suggest  that  the  linear  growth  at  resonance  will 
be  terminated  by  either  dissipative  effects  or  decay  of  the  internal  wave.  How¬ 
ever,  analysis  in  the  following  sections  indicate  that  wave  reflections  may 
also  determine  the  effective  interaction  time. 

The  mathematical  reason  for  the  linear  growth  at  resonance  in  equ.  (2.28) 
arises  basically  from  the  introduction  of  the  Born  approximation  in  the  calcu¬ 
lation  of  the  right  hand  side  of  equ.  (2.14)»  where  the  non-linear  terms  were 
replaced  by  products  of  the  unperturbed  surface  and  internal  waves,  i.e., 
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As  the  Interaction  proceeds  the  surface  wave  may  be  altered,  result¬ 


ing  in  failure  of  the  Born  approximation.  It  is  argued  later  in  this  paper 
that  this  is  indeed  what  does  occur  and  that  the  basic  process  is  one  of  "col¬ 
lision  or  scattering  (elastic  or  inelastic)  of  a  surface  wave  from  the  surface 


current  pattern.  For  most  of  the  strongly  interacting  "collisions"  the  effec¬ 


tive  interaction  time  is  determined  by  the  kinematics  of  this  collision  or 


scattering  process.  In  order  to  treat  the  near  resonance  waves  it  is  necessary 
to  proceed  beyond  the  Bom  approximation  and  follow  the  alteration  of  the 

surface  wave  as  it  interacts  with  the  current  pattern.  This  is  done  in  the 
following  section. 
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3.  WAVESCATTER  AND  REFLECTION 


We  start,  as  did  Rosenbluth  (1971)  and  Longuet-Higgins/Stewart  (1960), 
from  Equations  (2.11)  and  (2.12);  that  is,  linear  equations  keeping  first  order 
interaction  terms  involving  a  prespecified,  non-uniform,  steady  surface  cur¬ 
rent.  Rosenbluth  obtains  full  vave  solutions  for  a  surface  wave  moving  in 
a  region  of  steady  sinusoidal  varying  current  and  shows  that  surface  waves 
can  be  trapped  over  the  troughs  of  the  internal  waves.  The  solutions  ex¬ 
hibit  close  similarity  to  those  of  a  quantum  mechanical  harmonic  oscillator 
and  there  is  an  analogy  between  surface  waves  trapped  in  regions  of  minimum 
surface  current  and  quantum  mechanical  particles  trapped  in  regions  of 
minimum  potential  energy.  We  will  treat  the  problem  of  the  impingement  and 
reflection  of  a  free  surface  wave  from  a  region  of  linearly  increasing  sur¬ 
face  current  and  show  that,  for  a  certain  range  of  wavelengths,  the  surface 

waves  will  be  reflected  by  a  travelling  current  pattern. 

3A.  Dispersion  Relations  for  a  steady  surface  current 

Consider  a  finite  region  of  surface  current  translating  with  uniform 
velocity  and  described  by  a  function  U(x-ct).  We  transform  to  a  coordinate 
system  moving  with  this  pattern  (£  =  x-ct)  and  thus  obtain 

3  _  3  3  3  '  3  ,, 

3t  3t  "  °  3C  ’  3x  3 £  (3,1) 

so  that  in  this  new  coordinate  system  the  surface  waves  in  Equation  (2.11) 
obey  the  equations 

[♦«  -  2c*5t  +  c\{  +  2U  (c)kt  -  ■*«} +  g\Lrface  ■  0  <3-2) 

and 

<pFr  +  <t>  =  0  .  (3.3) 

££  yy 
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Taking  y  to  be  positive  upwards,  we  may  write  the  general  solution  of  Eq . (3.2) 
as  a  superposition  of  two  waves  [for  surface  wavelengths  «  depth  of  fluid  CD)] 


4>  = 

♦<R>  +  ♦«•> 

(3.4) 

where 

,(R) 

f  ky  ikx  a 

Je  ak(t)  e  dk 

(3.5a) 

and 

*a)  * 

f eky  ak(t)  e”ikX  dk  , 

(3.5b) 

which  are  similar 

in  structure 

to  equ.  (2.19).  The  quantities 

represent 

waves  traveling  to  the  right  and  waves  to  the  left  in  this  moving 

coordinate  system.  According  to  Eq.  (3.5),  we  may  write 


,«>  . 
y  x 

(3.6) 

,<L>  -  i*(L>  . 

y  x 

(3.7) 

It  is  to  be  noted  that  "traveling  to  the  right"  and  "traveling  to  the  left" 
refers  to  the  phase  velocity  of  the  surface  wave.  It  was  pointed  out 
in  section  2  that  for  this  surface-internal  wave  problem,  the  most  interesting 
phenomena  occur  when  the  group  velocity  of  the  surface  wave  is  closely  matched 
to  the  phase  velocity  of  the  internal  wave  equ.  (2.26)  .  Since  the  surface  waves 
are  dispersive,  it  is  possible  (and  this  will- be  the  condition  of  most  interest 
to  us)  to  have  a  situation  where  the  phase  velocities  are  in  the  same,  but  the 
group  velocities  are  in  opposite,  directions  for  two  coupled  waves.  It  occurs 
that  for  a  given  frequency  there  are  four  coupled  waves  that  appear  in  this  moving 
coordinate  system.  Two  of  these  waves  will  be  shown  to  have  phase  velocities  to 
the  right,  in  the  direction  of  the  iaternal  wave  (R  waves),  and  two  have  phase 
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velocities  in  the  opposite  direction  (L  waves).  For  the  two  R  waves,  one  turns 
out  to  have  a  group  velocity  to  the  right  (which  we  refer  to  as  a  (+)  wave: 

[,>(+K)l,  and  the  other  has  a  negative  group  velocity  [a  (-)  wave:  *  ] .  The 

two  L  waves  will  similarly  be  noted  as  (+)  and  (-)  waves  although  for  tl  esc 

the  group  velocities  have  the  same  sign. 

To  demonstrate  the  nature  of  these  waves,  we  consider  an  internal  wave 

packet  of  finite  size  such  as  the  one  shown  in  Fig.  3a.  A  packet  of  surface 
gravity  waves  propagating  slightly  slower  (group  velocity)  than  the  iirternrl 
wave  will  be  overtaken  and,  in  the  frame  moving  with  the  internal  wave,  these 
wave  packets  appear  to  impinge  from  the  right,  whereas  faster  surface  waves  impinge 
from  the  left.  We  shall  examine  the  wave  field  for  a  single  monochromatic 
surface  wave  incident  on  the  surface  current  region  (the  internal  wave  packet) 
from  the  left  [see  Fig.  3b]  i.e.,  one  moving  faster  in  the  moving  frame  of  refer¬ 
ence  than  the  internal  wave.  Since  in  the  moving  frame  the  surface  current 
field  is  steady,  the  time  dependence  of  the  surface  wave  potential  is  of  the 
form  e”iu>t  and  Eq.  (3.2)  becomes 

-u2*  +  i2cw*  +  c24>^  +  2U(0[-iw$5  -  +  g$y  =  0  •  (3,8) 

L  R 

This  oscillatory  time  dependence  occurs  in  the  coefficient  a^t)  and  a^(t)  in 
Eqs .  (3.5a)  and  (3.5b)  . 

Before  treating  the  case  of  a  space  variable  surface  current,  we  examine 
the  dispersion  relation  when  the  current  is  uniform,  i.e.,  U(£)  = Uq  [constant].  To 
study  this  case  we  will  reduce  Eq.  (3.8)  to  dimensionless  form  by  introducing  a 
characteristic  wave  number  (k*)  and  frequency  (u>*),  i.e.,  the  values  for  a  sur¬ 
face  wave  traveling  at  the  same  velocity  (phase)  as  the  internal  wave;  so  that 

U  ( 3  9) 

n  -  uj/w*  ;  K  m  k/k*  ;  8  =  o/c 
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Figure  3a,  Wave  Kinematics  Before  Interaction 


Figure  3b,  Wave  Kinematics  After  Interaction 
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where  a*  =  g/c  and  k*  =  g/c2.  Using  the  relation  (9/3?  -  ik)  from  Equations 
(3.1)  and  (3.6)  we  can  write  the  dispersion  relations  for  Equation  (3.8)  as 

-n2  -  2n*  -  k2  +  26 |n  +  +  <  =  0  for  the  (R>  waves  (3.10a) 

and 

_n2  +  2hk  -  K2  +  26  j-n  +  k]k  +  <  *  0  for  the  (L)  waves  (3.10b) 

thereby  giving  the  four  waves  discussed  earlier.  We  need  only  follow  the 
formalism  for  the  (R)  waves  since,  according  to  Equation  (3.10) 

4> (L)  (n)  =  <fr(R)(-h) 

and  for  these  waves 

1  -  2nd  -  8)  +  -  An (1  -  8)  + 

i  = _ . _ I _ •  _ -  ;  n  >  0  .  (3.11a) 

2(1  -  26) 


To  terms  linear  in  6  this  expression  agrees  with  the  one  obtained  more 
Intuitively  by  replacing  the  frequency  u>  in  the  dispersion  relation  calcu¬ 
lated  in  the  ocean  frame  (u>  =  ±/gk  )  by  the  value  expected  if  the  waves  are 
being  convected  at  a  velocity  (u-c),  i.e.,  u>  -*■  <*>  +  (c-u)k.  In  this  case 
the  dispersion  relation  would  become 

a)  +  (c  -  u)k  =  ±  \[gk 

or 

n  +  (1  -  6)<  K  ±  V* 

Solving  for  k  yields 
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K 


(3.11b) 


1  -  2n (l  -  6)  ±  -  And  -  6) 

2(1  -  6) 2 

which  agrees  with  Equation  (3.11a)  to  terms  linear  In  6.  Since  the  present 
analysis  applies  only  for  values  of  6  small  compared  to  unity,  we  will 
use  (3.11a)  or  (3.11b)  interchangeably. 


In  Figure  4  we  have  plotted  the  dispersion  relation  in  terms  of  the 
reduced  variables  n  and  <  .  We  are  particularly  interested  in  regions 
near  the  turning  point  where  the  group  velocity  vanishes  and  the  surface 
wave  nacket  travels  with  a  velocity  near  that  of  the  internal  wave  phase. 

In  this  region  the  two  R  waves  are  strongly  coupled.  However,  as  is  dis¬ 
cussed  later,  the  coupling  to  the  L  waves  is  weak  and  we  ignore  it  entirely 
here.  Thus,  we  delete  the  L,R  superscripts  and  refer  hereafter  only  to 
(+)  and  (-)  waves,  (both  of  which  will  always  be  R  waves). 


We  now  return  to  the  problem  of  a  finite  region  of  surface  current. 

We  are  treating  the  case  of  an  overtaking  packet  incident  on  this  region 
from  the  left.  This  wave  [a  (+)  wave],  in  general,  will  be  partially  trans 
mitted  (as  a  (+)  wave]  by  the  current  region  and  partially  reflected  [as  a 
(-)  wave].  if  we  define  r  as  the  coefficient  of  the  reflected  wave  and  t 

as  the  coefficient  of  the  transmitted  wave,  then  we  may  write  the  potential 
far  to  the  left  of  this  interaction  region  as 


i 

r  e 


(3.12a) 


The  meaning  of  reflected  and  transmitted  will  become  clear  later  as  we 
evelop  the  concept  of  the  surface  wave  packet  later  in  this  section. 
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where  $o  is  the  amplitude  of  the  incident  wave.  Further,  there  will  be 
only  a  transmitted  wave  far  to  the  right  of  the  region: 


t  e 


(3.12b) 


for  £ -*■  +<».  The  wavenumbers  and  ^  are  defined  in  the  undisturbed 

00  OD 

region  of  the  surface,  i.e.,  where  6=0.  One  may  therefore  use  the  dis¬ 
persion  relations  from  Equation  (3.10)  to  obtain  for  the  R  waves 


-  f  (*  "  Vl  -  *n  ) 

or 

2oic  -  4  uic 

g  ▼  ”  g 


(3.13a) 


(3.13b) 


If  ai  >g/4c,  therefore,  there  are  no  right  propagating  waves  in  the  frame  of  the 

moving  internal  wave.  The  distortion  or  reflection  of  the  surface  wave  as  it 

travels  through  the  current  region  depends  critically  on  the  variation  of  the 

(+) 

wave  number  k  with  changes  in  U  . 

It  is  useful  to  attempt  to  interpret  the  concepts  of  Ball  and  of 
Longuet-Higglnc /Stewart  in  terms  of  the  scattering  or  reflection  mechanism.  In 
the  reference  fvame  moving  with  the  internal  wave  the  strongly  interacting 
(reflecting)  free  waves  form  a  standing  wave  pattern  on  one  side  of  the  reflection 
point  (see  Fig.  5).  This  standing  wave  consists  of  two  wave  trains  having  (in 
this  frame)  equal  but  opposite  group  velocities  but  phase  velocities  in  the 
same  direction.  Although  the  mean  group  velocity  is  zero  (so  that  the  standing 
wave  pattern  does  not  move)  the  mean  phase  velocity  is  non-zero.  Individual 
wavelets  (or  points  of  constant  phase)  appear  to  move  through  the  standing  wave 
pattern.  For  example,  for  a  left  incident  wave,  a  wavelet  will  appear  to  be 
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Figure  5a.  Currents  Above  the  Internal  Wave 


Figure  5b.  Untrapped  Wave  Reflecting  from  the  Crest  Region 
of  the  Internal  Wave 


Figure  5c.  Waves  Trapped  in  the  Trough  of  an  Internal  Wave 


horn  on  the  left  side  of  the  standing  wave  region,  to  travel  to  the  right 
through  the  pattern,  its  amplitude  undergoing  a  quasi-periodic  modulation 
by  the  standing  wave  and  finally  dying  or  disappearing  at  the  front  edge  of 
the  standing  wave  region,  i.e.,  the  reflection  point.  According  to  the  Born 
approximation  of  Ball  and  Hartle/Zachariason ,  the  stationary  current  pat¬ 
tern  induces  a  secondary  wave  that,  together  with  the  incident  wave,  cre¬ 
ates  a  stationary  standing  wave  pattern.  This  induced  wave  must  therefore 
have  a  group  velocity  equal  and  opposite  to  the  incident  wave  (since  the 
mean  group  velocity  is  equal  to  the  phase  velocity  of  the  standing  wave, 
i.e.,  to  zero).  Thus,  it  should  be  possible  to  interpret  the  Born  approx¬ 
imation  results  near  resonance  in  terms  of  the  creation  of  a  scattered 
wave  whose  wavenumber  and  frequency  are  dependent  on  the  difference  between 
those  of  the  incident  wave  and  that  of  a  perfectly  resonant  wave  (i.e., 

one  for  which  c  *  c).* 
g 

3B;  Scattering  at  a  Boundary 

The  preceeding  section  was  concerned  with  the  long  range  effects  of  the 
internal  current  on  the  surface  waves.  This  section  will  study  the  surface 
waves  in  the  interaction  region.  To  this  end,  we  note  that  when  the  inter¬ 
nal  current  increases  over  a  limited  region  of  space  two  effects  contribute 
to  the  observables  of  interest:  1)  the  variation  of  the  wave  slope  as  it 
passes  through  the  surface  current  region,  and  2)  the  generation  of  reflec¬ 
ted  waves  at  regions  where  the  wavenumber  changes  rapidly. 

*  Although  a  mechanism  allowing  for  depletion  of  the  incident  wave  would 
probably  be  required  to  get  reasonable  results  at  later  times. 
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In  order  to  solve  this  latter  problem  for  somewhat  arbitrary  u(«,  we 
modify  the  formulation  of  the  problem  slightly.  We  divide  the  .urface 
current  region  into  a  number  of  sections,  each  having  a  width  small  com¬ 
pared  to  the  wavelength  of  the  internal  wave  (Fig.  6).  We  index  the  posi¬ 
tions  along  the  internal  wave  by  5j  and  let  the  interval  between  any  two 
steps  1  and  0+1)  be  A5j.  An  incident  wave  will  generate  both  transmitted 


and  reflected  waves  as  it  impinges  on  each  of  these  steps.  Within  Mch  .action 
the  .urf.ce  current  is  taken  to  be  uniform  .„  thet  the  w.v.  potential  in  the  re¬ 


gion  between  Ej  and  fj+J  has  the  general  form 


♦j  -  Aj  exp[l(kW?-„t)J  +  Bj  «wp[i(kWE-„t)j 


(3.14) 


here  kj  and  are  determined  from  the  dispersion  relation  for  the  JtI' 

slab.  The  values  of  Aj  and  Bj  are  found  by  matching  boundary  conditions  at 
all  the  interfaces  between  the  Individual  sections  (see  Figure  7).  The  num- 
ber  of  independent  component  waves  that  need  to  be  allowed  for  in  Equation 
(3.14)  is  two  in  „„r  present  case  (gravity  waves).  In  general  this  number 
is  equal  to  the  number  of  allowable  propagating  wavenumbers  that  have  the 
same  frequency.  In  the  capillary-gravity  regime  where  surface  tension  ef¬ 
fects  are  important,  it  is  possible  to  have  these  intersections  of  the  dis¬ 
persion  curve  for  a  given  frequency  „  (see  Figure  4),  thus  requiring  three 
independent  waves  in  Equation  (3.14).  We  carry  out  the  present  formalism 

for  gravity  waves  only.  The  generalisation  to  the  capillary-gravity  regime 
will  be  made  subsequently. 
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Figure  6,  Finite  Dif 
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Fig.  7a.  LEFT  INCIDENT  WAVE. 
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Fig.  7b.  RIGHT  INCIDENT  WAVE. 
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The  two  cases  to  be  considered  are  ahown  in  Figure  7:  an  incident  wave 
traveling  from  left  to  right  (+)  of  unit  amplitude  is  shown  in  Figure  (la) 
while  a  wave  of  unit  amplitude  incident  from  right  (-)  is  shown  in  Figure 
(7b).  Note  that  we  are  discussing  the  R-solutions  of  Equation  (3.10a), 
which  has  been  decomposed  into  these  (+)  and  (-)  waves  in  the  moving  coor¬ 
dinate  system  of  the  internal  wave.  Figure  7  is  actually  a  blow-up  of  the 
step  regions  shown  in  Figure  6.  Since  we  are  treating  quantities  discretely 

in  this  section,  both  the  reflection  (r  )  and  transmission  (t  )  coefficients 

J  J 

for  the  waves  carry  an  index  to  distinguish  the  boundary  to  which  they  refer, 
i.e.,  point  in  space.  An  additional  index  of  (+)  and  (— )  is  carried  by 
these  coefficients  as  well  as  the  wavenumbers  to  distinguish  the  direction 
(group  veJocity)  of  the  incident  wave  with  which  they  are  associated. 


The  continuity  of  the  amplitude  and  slope  of  the  waves  at  the  j1"*1 
boundary,  between  regions  A  and  B  in  Figure  7  requires  that 


4>i(C  =  Cj  )  +  ^  =  Cj  )  ” 


and 


♦  t  dC 


1  d*t 

♦_  dC 

♦  *  dC 

«rn 

I 

r 

(3.15) 


(3.16) 


where  the  subscripts  refer  to  the  incident  (i),  reflected  (r)  and  transmitted 
(t)  waves.  For  an  incident  wave  traveling  from  left  to  right  (Figure  7a), 
we  have  from  Equations  (3.15)  and  (3.16)  the  continuity  conditions  that 


36 


(3.17) 


and 


1  +  r 


(+) 


+  r 


(+) 


(3.18) 


Similarly,  for  an  incident  wave  traveling  from  right  to  left  [Fig.  7b],  we  obtain 


1 


+  r 


(-) 


(3.19) 


(-)+rf-)k^+)  =  tf")  kf")  .  (3-20> 

B  J  B  J  A 


The  reflection  and  transmission  amplitudes  across  the  boundary  may 
be  evaluated  in  terms  of  the  wave  numbers  in  the  two  regions  by  using  Equations 
(3.17)  -  (3.20), 


tf+)  -  1 

J  1 

'k(+) 

A 

-  »n/ 

-ri  • 

(3.21a) 

[k(+) 

[a 

-  <+)\/ 

k+) 

-ri  • 

(3.21b) 

J 

k('} 

kB 

-  kB+>!/ 

k‘i 

-n  • 

(3.22a) 

and 

r<'>  = 

J 

k<-> 
l  B 

-  ^i/ 

-■41  ■ 

(3.22b) 

The  value  of  obtaining  these  expressions  will  become  evident  when  we  allow  the 
interval  between  steps  to  become  very  small.  It  will  then  be  possible  to  represent 

these  differences  in  wavenumbers  in  terms  of  derivatives.  This  construction  is 
carried  out  later  in  this  section  and  is  also  used  in  Section  4. 
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T°  br0“de"  °Ur  ^  a  wave  of  unit  amplitude 

inoidaat  fro.  the  left  on  a  current  pattern  of  finite  extent  [see  Figure  8,. 

»  the  superposition  of  the  indent  -a  reflected  wave  to  the  left  of  the  re- 
transmitted  wave  far  to  the  right  [see  Equations  (3.12a)  and  (3.12b)] 

To  formalize  these  statements  we  define  a  two  vector  *  whos 

°  vector  y,  whose  components  are  the 

local  values  of  the  (+)  and  (-)  waves 


<P 


We  also  define  a 
wave  vector  (-«) 
right  hand  side 


scattering  matrix 
lar  to  the  left  of 


S  with  elements  a  which  relates  the 
the  current  region  to  the  vector  on  the 


O  If) 


In  order  to  evaluate  the  elements  of  S 

»  we  consider  two  hypothetical  scat¬ 
tering  experiments.  Suppose  firsr  a  rj.\ 

PPose  firSt  8  wave  unit  amplitude  is  incident  on 

the  r.glon  from  the  left  (Figure  s)_  ^  ^iitude  ^  _  is 

that  of  the  reflected  wave  [which  is  a  (-)  wave,  i,  de£lned  as  R(+, , 
Than  we  must  have 


0  J 


(3.23) 
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R<+)  =  _s  /  s 

K  S21  22 


(3.24a) 


T  “  S11  "  S12  S21/S22 


(3.24b) 


In  a  similar  manner,  for  a  wave  of  unit  amplitude  incident  from  the  right 
(Figure  8)  with  reflected  and  transmitted  amplitudes  ^  and  T  ,  respec¬ 


tively,  we  may  write. 


’Sll  S12"  ° 

*_S21  S22.  -T  - 


(3.25) 


S12/S22 


(3.26a) 


=  1/s 


22  • 


(3.26b) 


The  matrix  S  is  the  scattering  matrix  for  the  entire  current  region. 

The  reflection  amplitudes  ^  and  R^  \  as  well  as  the  transmitted  wav° 
amplitudes  T(+)  and  T*-)  also  refer  the.  entire  surface  current  region.  We 
may  use  Equations  (3.24)  and  (3.26)  to  express  the  elements  of  the  scattering 
matrix  in  terms  of  these  "measured"  reflection  and  transmission  coefficients: 


S22  "  1/T  • 
s, 2  =  R(")/T(";  * 

s21  =  -R(+)/T(_)  , 

su  -  T(+)  -  R<’>  R(+)/T<-> 


(3.27a) 

(3.27b) 


(3.27c) 


(3. 27d) 


These  elements  may  now  be  applied  to  the  analysis  of  the  characteristic  of 
the  jth  boundary  which  initiated  this  subsection. 
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Wo  <IH  l,u*  •'  two-component  vector  *  made  up  of  the  amplitudes  of  the 
right  [<^  ^  ]  and  Left  ^  ]  travelling  waves  at  a  point  just  to  the  right 


of  the  j^’  interface: 


r  *  w 
.(-) 


(3.28) 


Using  a  propagation  matrix  Pj  and  a  reflection  or  scattering  matrix  Sj  ,  we 
may  write  the  vector  at  the  (j+l)th  step  in  terms  of  its  value  at  the  jth 
step,  by 


Vi  '  pj  '  sj  '  *j 

,  th 


(3.29) 


where  the  waves  incident  on  thfe  j  boundary  are  scattered  by  the  matrix 
Sj  and  propagated  across  the  unLform  interval  from  to  by  the  propa¬ 

gation  matrix  .  If  the  entire  current  region  has  been  divided  into  N 
steps,  then  the  wave  amplitude  ^  to  the  right  of  this  region  is  related 
to  the  amplitude  to  the  left  of  the  region,  according  to 


ip 


or 


N+l 


\’+l 


fn  v 

j=i  f 


S  0 


1  ’ 


where  S  is  now  the  scattering  matrix  for  the  entire  region: 

N 

s  -  n  ?!  •  sj 

j=i 


(3.30) 


(3.31) 


(3.32) 
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Thus  the  matrix  S  introduced  in  Equation  (3.11)  for  a  wave  crossing  an  arbi¬ 
trary  region  can  b.  thought  of  as  being  constructed  from  a  sequence  of  weak 
scatterings. 

A  little  thought  will  convince  one  that  the  reflection  (scattering)  matrix 

S.1  3t  Cach  °f  the  steps  has  the  s*ne  structure  as  the  scattering  matrix  S  for 
the  entire  region.  We  may  therefore,  write  the  matrix  Sj  in  terms  of  the 
reflection  and  transmission  coefficients  rj  and  tj  as  follows: 


ft  (+>  .  /+>»,.(-> 


• '  '  r  '  # /t ' 

J  j  /Cj 


r<->/t(-) 

J  /CJ 


i/t 


(-) 


(3.33) 


The  propagation  matrix  Pj  carries  the  translation  information  for  the  process 
and  may  be  written  as 


If  we  assume  that  the  reflections  described  by  Equation  (3.31)  are  weak, 
then  it  is  possible  to  write  the  reflection  matrix  above  [Equation  (3.33)] 
in  an  approximate  form.  For  weakly  reflected  waves  we  can  neglect  all  terms 
in  Sj  which  are  .econd  order  in  the  reflection  coefficients.  r<+)  and  r<">. 
Also,  we  may  expand  the  inverse  of  the  (-)  transmission  coefficient  as, 


1 


1  +  r 


(-) 

j 


(3.35) 


) 


> 
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Using  Equation  (3.35)  and  neglecting  terms  quadratic  in  the  reflection  coef¬ 
ficients,  we  write  the  scattering  matrix  as, 


where 


r  (+) 

(-)l 

r  J 

rj 

(+) 

(-) 

l~ri 

u 

1 

(3.36) 


(3.37) 


2 

where  is  considered  negligible  in  this  approximation.  This  Implies  that 
there  is  only  a  single  scattering  from  each  of  the  steps  and  that  the  reflected 
surface  wave  is  a  linear  superposition  of  each  of  these  small  reflections. 


Let  us  consider  the  generation  of  these  reflected  waves  by  Equation 
(3.32).  If  we  impose  our  weak  reflection  condition,  Equation  (3.36),  then  all 
terms  of  order  and  higher  in  this  product  may  be  neglected.  A  general  term 
which  would  contribute  to  the  scattering  matrix  would,  therefore,  be  linear 
in  R j ,  i.e. ,  be  of  the  form 


3  p  p  d  p 

N  N-l  •  •  •  j  "j  rj-] 


•  •  *  P2  P1 


(3.38) 


Physically,  Equation  (3,38)  means  that  the  Incident  wave  will  propagate  freely 
through  the  first  J-l  steps  of  the  steady  current  region;  at  the  jth  step  it 
will  interact  weakly  with  the  internal  wave,  as  specified  by  the  matrix  Rj, 
and  then  continue  to  propagate  without  further  interaction. 


One  can  easily  show  that  the  product  of  propagation  matrices  preceding 
the  interacting  in  Equation  (3,38)  can  be  written  in  the  continuous  limit  as, 
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p  p  .  .  .  p  p 
J-l  J-2  2  1 


expji /'  k(+)(C')dC' 


exp  k(_)(C')dC' 


(3.39) 


by  using  the  definition  of  the  P's  [Equation  (3.34)].  It  is  also  clear  that 
in  general,  the  amplitude  of  the  incident  wave  is  being  modulated  by  the  in¬ 
ternal  wave  packet. 

3C.  Differential  Formulation  of  the  Propagation  Equations 


In  this  section  we  will  construct  a  set  of  differential  equations  to 
govern  the  variations  of  the  forward  and  backward  propagating  waves  [4>^ 
and  by  taking  the  limit  of  very  small  S  intervals  in  the  propagation 

equations  [Equations  (3.21)  to  (3.22)].  In  the  limit  of  small  “  A£, 

the  transmission  and  reflection  amplitudes  become, 


(3.40) 


,<+) 

AS 

dK<-> 

j 

[k<+>  -  K<-)| 

dS 

.(-) 

4{ 

d<(-> 

j 

|«<+>  -  K<— >] 

d? 

(-) 

j 

=  1 

+  r<-> 

J 

.<+) 

=  1 

+  r<+>  . 

(3.41) 


(3.42) 


tj(+)  =  1  +  r<+)  .  (3.43) 

Equations  (3.40)  and  (3.41)  will  be  derived  in  Section  4,  where  a  complete 
discussion  of  their  characteristics  will  be  given. 
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The  dispersion  relation  [Equation  (3.11)]  (neglecting  quadratic,  terms  in 
0)  may  be  written  in  the  form 


where 


K±  =  K  +  K 

-<■  2 


k  =  l  -  2nd  -  01 
2(1-20) 


(3.44) 


(3.45a) 


v l  -  4n(l  -  8) 
2(1  -20) 


(3.45b) 


For  very  small  values  of  the  difference  5j+1  -  ^  the  scattering  matrix 

for  weak  reflections  [Equation  (3.36)],  using  Equations  (3.40)  to  (3.45),  can 


be  rewritten  as 


-  f1  °]  +  1  xl+  f’1  1 

[o  ij  _Xj  d5  2k2  1 


(3.46) 


Similarly,  the  propagation  matrix  given  by  Equation  (3.34)  may  be  expanded 


to  obtain 


,  .  f 1  °  1  r^Oei-KjWtj 

j  + 

1°  1  J  0  ik*(K1+  K2)lli 


From  the  recursion  relation  [Equation  (3.29)]  we  have 

d*1  /  X 

ar  *  Vi  -  *j  ■  (pjSj  -  i)*j 


(3.47) 


(3.48) 


where  ^  i8  defined  as  the  two  component  wave  given  by  Equation  (3.28).  In¬ 
troducing  Equations  (3.46)  and  (3.47)  into  Equation  (3.48)  and  keeping  only 
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terms  which  are  linear  in  ,  we  obtain  the  following  coupled  equations  for 
the  wave  components  ♦  ^  and 

5T—  ‘  +  S73rl*<+>+*<‘>l  -S-1T  [♦<+)‘*<'))  <3-W 


and 


50) 


in  the  continuous  limit.  These  equations  may  be  simplified  if  we  introduce 
the  sum  and  difference  amplitudes  S  and  D: 


and 


S  -  |$<+>+  expj-iy*k*  dgj 

D  -  exp|-iy"k*  d$j 


(3.51) 


(3.52) 


By  adding  and  subtracting  Equations  (3.49)  and  (3.50)  and  using  Equations 
(3.51)  and  (3.52)  we  obtain, 


dS 

dt 


and 


-ik*  <2  D 


<!°  .  -_L±iD  +  [.L±I_ 

d?  k2  de  D*  k2  d? 


ik*  k. 


(3.53) 


(3.54) 


We  may  eliminate  D  from  these  latter  equations  with  a  little  algebraic  manip¬ 
ulation  and  obtain 


2  d*l 

(k*»c2)z  +  ik* 


0 


(3.55) 


as  the  equation  for  the  wave  envelope. 
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3D.  Solutions  In  the  Neighborhood  of  a  Turning  Point 

Near  <1  -  0  ,  Equation  (3.49)  and  (3.50)  are  singular.  Near  these 
singular  regions  (turning  points),  we  may  approximate  the  current  pattern 
by  a  linear  form.  Equation  (3.55),  familiar  to  those  concerned  with  electro¬ 
magnetic  propagation  in  a  medium  having  a  variable  refractive  index,  has  an 
explicit  solution  when  the  current  (■  c0  ■  -ce)  varies  linearly  with 
position.  To  see  this,  let  us  define  a  critical  value  for  the  adverse 
current  e_  ,  i.e. ,  e0  ■  -i - 1  ,  and  define 


e  **  cq(1+5/A)  for  -A  <  £  <  0  . 


(3.56) 


Using  the  definitions  of  and  from  Equation  (3.45)  in  conjunction  with 

Equation  (3.56),  we  find  that  the  imaginary  term  in  Equation  (3.55)  is 

2 

negligible  in  comparison  to  the  (k*^)  term  when 


Id  »  i-  -  *1  . 

k*  * 


(3.5’) 


With  this  condition  Equation  (3.55)  reduces  to 


(3.58) 


which  is  valid  when  the  wave  envelope  is  much  more  than  a  characteristic  wave¬ 
length  away  from  the  reflection  point. 


Equation  (3.58)  is  the  equation  for  the  Airy  integral  in  terms  of  the 


variable  ( : 


2  \l/3 


<  -  (ti2) 


(3.59) 


that  is. 


-  CS  -  0  . 


(3.60) 


The  form  of  the  solution  is  that  of  a  varying  amplitude  standing  wave. 


The  wave  amplitudes  which  are  solutions  of  Equations  (3.49)  and  (3.50) 
are  expressed  in  terms  of  the  Airy  integral  [Ai(c)l  according  to 


\ 

where 


(3.61) 

(3.62) 

(3.63) 


from  Equation  (3.60),  and 


D 


dC 


(3.64) 


from  Equation  (3.53).  Here,  £  is  expressed  in  terms  of  £  in  Equation  (3.59) 
and  is  a  constant.  The  constant  can  be  determined  by  examining  the 
asymptotic  form  of  Ai(c),  which  for  C  <<  0  is, 


Ai(c)  ~ 


(3.65) 


Thus,  using  Equation  (3.63),  we  obtain  for  the  asymptotic  form  of  the  total 
wave  envelope, 

S  '  A  |c|1/4  s1"(^!|3/‘  +  ?)  ■  <3'66) 

Using  Equation  (3.66)  in  conjunction  with  Equation  (3.64)  we  find 


(3.67) 


as  the  asymptotic  expression  for  the  wave  difference.  Normalized  to  a  unit 
incident  amplitude  [<f>^  =1]  at  £--A,  the  asymptotic  form  of  the  incident 
and  reflected  waves  are 
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(ff  ^[A  «-iV?ui3/2]-fj 
(ity)7  c*pjik* [./"ki de + t ui  ]+i^j 


The  velocity  potential  of  the  total  wave  has  the  modulated  envelope: 


as  does  the  wave  height  yg.  The  wave  height  may  be  obtained  from  the 


>,  ■/&'■•  /•,  * 

which  using  Equations  (3.5)  and  (3.6)  yields 


±L+2 

.<+> 


and 


For  an  incident  wave  whose  peak  amplitude  at  £  ■  -  A  is  AY,  the  wave 
as  a  function  of  position  is  the  real  part  of 


sin  z 


exp 


lk  *fKi 


d£  - 


where 


or  in  other  words 
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(3.68) 

(3.69) 

(3.70) 

equation 

(3.71) 

height 

iiotl 


(3.72) 


4.  approximate  solutions 


We  now  use  the  severing  Eorm.ll..  of  the  preceding  section  to  examine 
the  approximate  change  In  the  surface  wave  as  It  Interact,  with  the  Internal 
current.  Using  again  the  notion  that  If  the  surface  current  varies  rapidly 
with  position,  reflected  waves  can  be  generated.  If  we  assume  that  these 
reflection,  are  weak,  then  It  Is  possible  to  write  the  reflection  matrix  In 
the  approximate  form  given  by  equation  (3.37).  These  assumptions  will  be  made 
in  the  following  section,  where  th.  reflected  wave  Is  calculated  In  the  Born 
approximation.  Section  4B  will  take  Into  account  the  variability  of  the  re¬ 
flected  wave  amplitude  and  discus,  the  elkon.l  (WKB)  solution  to  the  problem. 

It  will  be  determined  that  the  WKB  solution,  agree  with  the  solutions  obtained 
for  the  coupled  differential  equations  [(3.49)  and  (3.50)1  In  Section  3C. 

4A.  The  Bom  Approximation 

We  will  now  continue  the  analysis  begun  in  Section  3B,  where  the  scatte 
ing  „f  a  wave  at  a  boundary  1.  discussed.  In  this  section  we  assume  that  the 
amplitude  of  the  Incident  wave  1.  constant  throughout  the  Interaction.  We 
will  further  assume  that  the  end  points  of  our  propagation  Interval  (^  and  ?„) 
fn  Equation  (3.30)  are  very  far  fro.  the  Interaction  region.  Incorporating  the 
single  reflection  hypothesis  In  Equation  (3.30).  taking  the  continuous  limit  of 
the  product  of  the  propagation  matrices  [Equation  (3.39)]  and  using  the  scat- 
taring  matrix  fro.  Motion  (3.36).  .«  may  «lt.  for  th.  r.flacted  wave. 


(U  exp  i 


-  k(')tt')|d5'! 


(4.1) 
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Equation  (4.1)  is  the  continuous  analogue  of  Equation  (3.38),  i.e.,  because 
there  is  only  a  single  scattering  from  each  of  the  steps  (Cj)*  the  reflected 
surface  wave  is  a  linear  superposition  of  each  of  these  small  reflections. 


We  may  write  Eq.  (4.1)  in  a  more  manageable  form  by  taking  the  limit  of  the 
series  as  an  integral  in  the  same  spirit  as  we  constructed  the  propagators  [Eq. 
(3.39)].  To  do  this,  consider  the  reflection  coefficient  in  Eq.  (3.21b),  where  the 
wavenumbers  and  k^  may  be  more  instructively  written  in  terms  of  the 

position  variable  £  .  We  delete  the  index  labeling  the  region  and  write 
the  difference  in  wavenumber  at  the  jt'1  step  as 


.(+) 


-5 


dk 


(+) 


■J+l  jl  dC 


(4.2) 


5  =  5, 


when  the  step  length  |Cj+1  -  Cj j  is  sufficiently  small.  Using  Eq.  (4.2) 
and  (3.21b),  we  may  write  the  reflection  coefficient  as  a  function  of  posi¬ 
tion 


rj(+)(Cj)  »'  - 


k(+)(Cj)-k<")(Cj) 


dk 

dC 


(+) 


5-C. 


(4.3) 


Further,  since  the  variation  in  wavenumber  with  position  is  dependent  on  the 
variation  in  the  surface  current  (“  8c)  with  position;  it  is  more  natural  to 
write  equ.  (4.3)  as 


r<+>» 


dk(+)dg 
dg  dC 


(4.4) 


By  substituting  Equation  (4.4)  into  Equation  (4.1)  we  may  express  the  reflected 
wave  in  terms  of  the  integral, 
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♦ref]  _  **P  ♦■■*/[«  W(e')  -  dc']  di<(+)  JB 

*o  '  x<+)(«)  -  K(->(e>  « 


d£ 


(4.5) 


where  the  continuous  limit  of  Equation  (4.1)  has  been  taken.  A  further  simpli¬ 
fication  can  be  made  near  the  resonance  condition  8  +  0,  n+1/4,  i.e.,  the  point 
where  <2  becomes  zero  in  Equation  (3.44): 

\1n_  r  exp{lk*/£dcj  ^ 

♦o  J  B7C  /rr^T 

where  k  =  [tc^  ^  ^]/2,  and  we  have  used  p2  —  1/16  in  an  expansion  of 

/l  -  4n  ■ 


In  the  Born  approximation  the  wavenumber  k  is  assumed  to  be  a  constant 


throughout  the  interaction  interval  and  may,  therefore,  be  removed  from  the 
integral  sign  as  follows 


ref  1 


8x  c  /l  -  4n 


i2k*tc£  dU 


d£ 


(4.6) 


We  may  integrate  Equation  (4.6)  by  parts  to  obtain, 


+ 


k* 

4/i  -  4n 


ei2k*<C 


d£  • 


(4.7) 


This  is  the  Bragg  condition  for  scattering,  where  the  reflection  from  the 
train  of  surface  currents  is  proportional  to  the  Fourier  component  of  the 
surface  current  distribution  at  the  "mean"  wavenumber  2k*i<  .  Since  at  the 
resonant  condition  the  surface  gravity  waves  are  typically  of  much  shorter 
wavelength  than  the  internal  waves,  this  contribution  to  the  reflection  is 
usually  very  small. 
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4B.  The  Eikonal  (WKB)  Approximation 


Of  much  greater  Importance  than  the  weak  Bragg  reflections  considered 
in  Section  4A  is  the  possibility  of  total  reflection  of  the  impinging  waves, 
and  associated  with  this,  the  possibility  of  trapped  waves.  To  analyze 
this  problem  we  will  use  the  WKB  approximation  in  the  slowly  varying  regions 
far  from  the  turning  points  and  match  theoe  asymptotic  forms  to  the  Airy 
integral  expressions  (Section  3D)  which  are  valid  in  the  vicinity  of  the 
turning  point. 

The  normal  WKB  approximation  (as  used  by  moat  workers)  permits  only  one 
wavenumber  to  characterize  a  wave  packet  and  thus  must  be  carefully  interpreted 
when  applied  to  the  problem  of  a  wave  packet  incident  on  a  region  having  reflecting 
properties.  In  the  coordinate  system  moving  with  the  pattern  of  surface  current 
produced  by  the  internal  wave,  there  are  four  different  propagating  waves  which 
may  coexist  at  the  same  frequency  (see  Fig.  2  ).  Two  of  these  can  be  strongly 
coupled  if  the  wave  frequency  is  close  to  the  resonant  value  (g/4c).  Thus,  when 
applying  the  WKB  approximation  to  this  problem  it  is  important  to  allow  for  at 
least  two  coexistent  waves.  The  effects  are  most  easily  deduced  from  the  full 
wave  equation. 

We  consider  first  the  distortion  of  the  wave  as  it  passes  through  a 
region  of  slowly  varying  surface  currents,  described  by  U (O •  Thus,  at  point  £, 
we  have  in  general  for  a  right  traveling  wave,  using  Eqs.  (3.1),  (3.6)  and  (3.8), 

-io24>  +  2ciu><(^  +  +  2U  (C)j-iu^  -  c$^J  -  ig<^  -  0.  (4.8) 

We  wish  to  use  the  WKB  approximation  to  describe  the  distortion  in  the 
interaction  region,  to  this  end  we  introduce  the  normalized  variables  from 
Section  2,  and  rewrite  Equation  (4.8)  as, 
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i  .  lk*  2n(i-g)  -1  a  _  n2-^  4-0 
4  ik  1  -  20  *5  i  -  20  ? 


(4.9) 


To  gain  insight  into  Equation  (4.9)  let  us  put  the  equation  in  the  fol¬ 


lowing  form, 


^ _  (P$)  +  GPj  ■  0 

dC2 


(4.10) 


Taking  the  prescribed  derivatives  and  equating  terms  to  those  ir  Equation  (4.9) 


yields  for  the  functions  P  and  G 


P  =  exp 


j-i /k*  <1  dC 


(4.11) 


G  -  (k*K,)  +  ik* 


(4.12) 


Therefore,  with  the  substitution  S-  P$,  we  obtain 


£*  +  [(k*«2)2  +  ik*  ^l]s  -  o 

nr  L 


(4.13) 


for  Equation  (4.9),  which  is  just  that  obtained  from  the  coupled  differential 
equations  in  Section  3C  [3.55],  when  P+  is  the  total  wave  envelope. 

In  the  eikonal  or  WKB  approximation  we  write  5  in  Equation  (4.13)  as, 


S  -  Ae 


(4.14) 


when  we  require  that  both  the  amplitude  (A)  and  phase  (T)  be  real.  Substituting 
Equation  (4.14)  into  (4.13)  and  equating  real  and  imaginary  parts  yields  the 


coupled  equations, 


<k*K2)2  -  r2*  -  0 


(4.15) 


ree*2re  VA  +  k*3T  "  °  * 


(4.16) 
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It  Is  convenient  to  define  ■  <  in  the  above  equations,  which  may  then  be 
written  as 


K 


2 


(k*K2)2 


(A. 17) 


and 


d 

dT 


(A2<) 


2 


(A. 18) 


When  <  is  slowly  varying  over  many  wavelengths,  then  we  may  neglect  the 
second  term  on  the  right  hand  side  of  Equation  (A. 17).  In  this  limit  (the 
elkonal  approximation) 


k  -  ±  k*  k2(5) 
and  ^ 

r  -  ±/kV,')  dc' 


(A. 19) 


A  second  approximation,  which  is  consistent  with  the  analytic  treatment  in 
Section  3C  is  to  neglect  the  contribution  of  d'^/'iC  to  Equation  (4.18).  In 
this  case  the  amplitude  varies  according  to 


AiiL. 

A(50) 


▼ 


) 


0_ 

CO 


which  results  from  the  approximation 
(A2k)  «  0  . 


(A. 20) 


Using  these  results  in  conjunction  with  Equations  (A.1A)  and  (4.11)  we  may 
write  the  envelope  of  the  wave  potential  in  the  form, 


iut 


(A. 21) 


where  is  the  initial  amplitude.  The  wave  height  at  a  point  €  with  reference 
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to  the  height  at  the  point  Cq  where  the  current  [U(C)1  vanishes,  is  given  by 


a/4 


JL  - 
y. 


exp 


ji J  Ic*[ic1(C#)  ±  x2(C')]dc'  - 


iut 


(4.22) 


where  and  <2  may  be  obtained  from  Equation  (3.45).  Equation  (4.21)  is  in 
essential  agreement  with  the  asymptotic  expression  in  Equation  (3.69)  when 


a  linear  current  profile  is  assumed  and  when  we  identify  the  +  sign  with  $ 


<♦) 


and  the  minus  sign  with  Thus  we  may  take  the  general  solution  to  be 

given  by 


.  .<♦>  +  .<-> 


where  $  ^  and  $  ^  are  given  by  Equation  (4.21),  except  in  the  vicinity  of 
the  turning  point.  In  this  latter  region,  we  must  use  the  full  wave  expression 
(Eqs.  3.63  and  3.65). 
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5.  DISCUSSION  AND  CONCLUSIONS 


This  section  is  divided  into  three  parts.  The  first  part  summarizes 
the  mathematical  structure  of  the  paper  stressing  the  various  assumptions 
employed.  A  descriptive  summary  of  the  physical  mechanisms  involved  in 
the  scattering  process  based  on  analogy  with  balls  rolling  on  slowly-varying 
potential  surfaces  is  provided  in  the  second  subsection.  The  third  sub¬ 
section  discusses  the  phenomenon  of  surface  modifications  produced  by  weak 
currents  and  presents  exemplary  calculations  of  the  modifications  expected 
for  a  'linear'  sea. 

5A.  Summary  of  the  Formal  Structure 

In  Figure  9  we  indicate  diagramatically  the  major  content  and  structure 
of  the  mathematical  analysis  developed  in  Sections  3  and  4.  The  general 
hydrodynamic  equation  for  the  ocean  surface 

+  (V*‘V*)  +  i  V**V(Vf.V#)  +  g#  1  -  0  (5.1a) 

t  2  ylsurface 

is  used  to  model  the  interaction  between  surface  gravity  waves  and  low  speed, 
low  amplitude  surface  currents.  This  is  done  by  assuming  the  velocity  potential 
to  be  of  the  form  +  where  U(x,t)  -  V*(o)  is  the  weak  surface  current  and 

$  is  the  velocity  potential  for  the  surface  gravity  waves.  Equation  (5.1)  is 
linearized  by  taking  the  ratio  of  the  gravity  wavelength  to  the  current  wave¬ 
length  (X  /A  )  <<  1  ,  so  that  in  one  horizontal  dimension  we  obtain 

♦  tt  +  S+y  +  2U(xit)4xt  *  0  .  (5.1b) 
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FIGURE  9,  Structure  of  the  Mathematical  Development 
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Equation  (5.1b)  has  a  first  order  interaction  term  involving  a  prespecified, 
non— uniform  surface  current. 

There  are  a  number  of  possible  causes  of  the  surface  current,  e.g., 
general  ocean  circulation,  swells,  river  inlets,  internal  waves,  etc.  If 
we  assume  that  the  current  is  the  result  of  an  internal  wave  propagating 
along  a  thermocline,  the  situation  is  unique  in  that,  in  the  coordinate 
system  travelling  with  the  internal  wave,  the  surf ace  current  pattern  is 
stationary  (so  long  as  the  internal  wave  is  non-dispersive) .  In  this 
coordinate  sytem  (C  -  x  -  ct)  there  exist  four  possible  co-existing  waves, 
two  of  which  have  a  positive  and  two  of  which  have  a  negative  phase  velocity 
[see  Figure  (4)].  The  two  waves  with  positive  phase  velocities  (right 
moving  or  R-waves),  have  a  dispersion  relation  (for  constant  current  UQ) 

-n2  -  2m  -  k2  +  2B0  +kJk  +  K  =  0  (5.2) 

in  terms  of  the  normalized  frequency  (n  -  <%>/g),  wavenumber  (tc  =  kc2/g) 
and  current  (S0  ■  bj/c).  These  two  waves  have  opposing  group  velocities 
(doj/dic  -  dn/dic  -  Cg) ,  however,  and  our  discussion  centers  on  the  inter¬ 
dependence  of  these  two  waves  and  waves).  In  the  case  of  a 

oteady  current  [ U ( O  ] ,  this  description  leads  in  a  natural  way  to  the  notion 
of  incident  and  reflected  waves  from  the  weak  current  region. 

From  the  above  dispersion  relation  [Equation  (5.2)]  it  is  possible  to  deter¬ 
mine  the  wave  frequency  in  terms  of  the  wavenumber  and  current  amplitude.  At  the 
point  where  the  group  velocity  of  the  R-wave  vanishes  (c  *  dn/d<  «  0)  the  surface 

O 

wave  travels  in  synchronism  with  the  internal  wave  and  the  and  $(_)  waves  are 

strongly  coupled.  This  condition  for  strong  coupling  can  be  expressed  in  terms  of 
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critical  values  for  the  current,  via  the  dispersion  relation.  These  concepts 
are  formulated  mathematically  in  Section  3B  by  considering  a  stepwise  constant 
current.  Each  step  is  treated  as  a  boundary  and  the  incident  wave  is  partially 
reflected  and  partially  transmitted  at  each  boundary.  By  requiring  continuity 
of  amplitude  and  slope  of  the  transmitted  and  reflected  waves  at  each  of  these 
steps,  a  scattering  matrix  for  the  current  region  may  be  developed.  The  con¬ 
tinuous  limit  of  these  equations  using  the  assumption  of  a  weakly  reflecting 
current  is  taken  in  Section  3C.  This  leads  to  a  set  of  coupled  partial  dif¬ 
ferential  equations  for  and  ^ .  The  coupling  is  dependent  on  both 

the  current  and  the  spatial  gradient  of  the  current. 

The  above  set  of  coupled  differential  equations  [Equations  (3.49)  and 
(3.50)]  may  be  reduced  to  a  single  equation  for  the  wave . envelope,  i.e., 

S  ~  ,  [Equation  (3.55)].  This  equation  becomes  an  Airy  integral 

for  a  linear  current  profile  when  the  distance  from  the  reflection  point  is 
much  greater  than  the  characteristic  wavelength  of  the  gravity  waves. 

In  Section  4  the  connection  between  the  scattering  approach  to  the  in¬ 
teraction  problem  is  compared  with  a  straightforward  treatment  of  the  linearized 
dynamic  equations  [Equation  (5.1b)].  A  little  algebra  establishes  that  in  the 
moving  coordinate  system 

.2Q 

+  (a  +io2)S  -  0  (5.3) 

d  r 
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where  S  is  the  wave  envelope  as  defined  in  Section  3.  The  scattering  theory, 
therefore,  yields  the  same  dynamic  equation  for  the  interaction  problem  as 
linearizing  Equation  (5.1). 

Equation  (5.3)  is  solved  in  the  WKB  solution  under  the  condition  that 
>>  oi^ ,  which  is  the  same  condition  used  in  Section  3.  The  WKB  solutions  are 
Identical  to  the  asymptotic  forms  obtained  in  Section  3  when  a  linear  profile 
for  the  current  is  assumed.  The  turning  points  in  the  WKB  solutions,  i.e., 
the  points  where  the  wave  amplitude  becomes  infinite,  are  the  spectral  or 
total  reflection  points  for  the  incident  wave.  This  is  discussed  more  fully 
in  Section  5C. 

5B.  Interpretation  in  terms  of  the  scattering  of  wave  packets. 

The  above  summary  of  the  formalism  has  shown  that  the  two  models  - 
(i)  the  linearized  hydrodynamic  equation  for  the  ocean  surface,  and  (ii)  the 
notion  of  weak  reflections  of  the  surface  gravity  waves  by  the  surface  current  ■ 
imply  the  same  dynamic  description  of  the  interactive  process.  The  importance 
of  the  present  description  of  the  process  lies  in  the  physical  insight  that 
it  provides.  We  see  that,  just  as  in  the  scattering  of  electromagnetic 
waves  from  a  region  of  varying  refractive  index  or  scattering  of  a  quantum 
mechanical  wave  (i.e.,  a  particle)  from  a  potential  well,  it  is  possible 
to  interpret  this  linearized  ocean  wave-current  interaction  model  in  terms 
of  scattering  theory. 
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Under  certain  circumstances  the  scattering  can  be  very  strong 
(resonant).  The  interaction  has  a  close  formal  analogy  with  scattering 
of  particles  by  a  shallow  potential  well.  In  order  to  give  a  physical 
description  of  the  type  of  resonance  involved,  we  describe  the  analogous 
resonance  which  occurs  in  scattering  of  particles  by  a  shallow  potential 
well  and  make  reference  to  the  text  where  appropriate. 

Think  of  a  tightly  stretched  rubber  sheet,  and  of  a  very  slight  hill 

travelling  at  constant  speed  (vb^j_)  across  this  sheet  (as  would  be  created  by 

a  moving  finger  pushing  the  sheet  up  just  slightly  from  below).  The  hill  is 

the  surface  current  distribution  (v  ■  c  .  )  created  by  the  internal  wave. 

hill  phase 

We  now  shoot  marbles  at  the  moving  hill.  The  marbleB  are  the  short  wavelength 
surface  gravity  waves  (vmarble  ■  c^) .  A  marble  moving  very  faBt  will  not  be 
much  affected  by  the  hill  except  for  feeling  it  as  a  brief  faint  bump.  A  marble 

moving  very  slowly  will  also  not  be  affected  very  much  since  the  hill  will 
overtake  it  and  pass  it  too  quickly  to  cauBe  much  motion.  However,  a  marble 
that  moves  almost  at  the  same  velocity  as  the  hill  feels  its  presence  for  a 
long  time  and  can  interact  strongly.  This  is  the  resonance.  For  low  hills  a 
perturbation  theory  could  be  used  to  describe  the  interaction,  but  the  results 
would  blow  up  when  v^^  ■  vmarble  (Equation  (2.26)],  giving  a  perturbation 
growing  linearly  in  time  [Equation  (2,28)]  (this  is  essentially  the  Hartle- 
Zachariason  calculation  for  the  wave-wave  interaction).  In  this  model 
(at  least  in  a  one-dimensional  model),  it  can  be  shown  that  there  if.  a  specific, 
sharply  bounded  range  of  marble  velocities  centered  around  v^^  for  which  the 
interaction  is  strong  and  that  outside  this  range  the  interaction  is  weak.  In 
this  velocity  range  the  marbles  lose  or  gain  energy  as  they  hit  or  are  hit  by 
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the  hill.  The  effect  is  aiost  easily  understood  by  viewing  it  in  a  coordinate 
system  moving  with  the  hill  (remember  this  is  one  dimension  so  that  the  hill  is 
more  like  a  long  ridge).  In  this  coordinate  system  the  hill  (internal  wave  cur¬ 
rent)  does  not  move.  Marbles  (short  gravities)  overtaking  the  hill  appear  now 
to  come  from  the  left  whereas  marbles  being  overtaken  by  the  moving  hill  appear 
to  impinge  from  the  right.  The  impact  velocity  is  ±  (vmr^e  -  v^^)  ■  ±|Av| 

(see  Section  1).  For  large  | Av |  the  marble  slows  down  slightly  as  it  climbs  the 
hill,  but  regains  its  original  velocity  as  it  rolls  down  the  other  side.  How¬ 
ever,  there  is  a  critical  value  of  Av(=  Av*)  ,  dependent  on  the  hill  height  (AQ) , 
below  which  marbles  cannot  get  over  the  hill  -  they  are  reflected.  In  the  ocean- 
fixed  coordinate  frame,  these  marbles  gain  or  lose  a  velocity  ±2Av.  The  interaction 
is  very  strong  for  marbles  having  velocities  in  this  range  (vhill  ±  Av*(Ao)), 
i.e.,  every  9uch  collision  results  in  a  (group)  velocity  change  of  ±2Av. 


In  particular,  it  is  to  be  noted  that  n£  marbles  which  have  a  velocity 
in  the  range  vhlll  ±  Av*(AQ)  are  ever  found  on  or  near  the  top  of  the  hill  - 
they  are  always  reflected  away. 

For  a  given  internal  wave  having  speed  c  and  amplitude  AQ,  there  is  a  specific 
range  (X^  to  X2)  wavelengths  of  surface  gravity  waves  which  are  strongly  re¬ 
flected  by  the  current,  pattern  of  the  Internal  wave.  That  is,  for  X^  <  X  <  *2» 
waves  obeying 


c 


<  Av*(AQ) 


cannot  exist  at  the  top  of  the  internal  wave.  Such  wavelengths  should  be  absent 
or  at  least  strongly  altered  in  the  wave  amplitude  power  spectrum  and  therefore 
also  in  the  radar  return  in  certain  regions  over  the  internal  wave. 
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5C.  Application  to  the  Mechanical  Generation  of  'Slicks' 

In  this  section  we  provide  estimates  of  the  range  of  frequencies  and 
wavenumbers  which  are  predicted  to  be  strongly  altered  by  a  given  internal 
wave.  Also,  examples  of  the  spatial  modulation  of  a  wave  spectrum  are 
evaluated  for  a  'linear'  sea. 

For  a  finite  packet  of  roughly  sinusoidal  currents  having  a  peak 
amplitude  A0,  the  frequencies  which  can  be  reflected,  i.e.,  for  which  there 
are  turning  points,  should  be  in  a  range 


that  is, 


u 


i 

4  c 


K  1 

JL  i 

4c  4  c2 


(5.5a) 


(5.5b) 


These  are  wave  frequencies  a.  seen  in  the  moving  coordinate  system.  The  cor¬ 
responding  wavenumber  range  (outside  the  current  region)  can  be  obtained  from 
the  dispersion  relation  for  the  unperturbed  ocean t 


kc‘ 


■  K  ■ 


L-Jl  ±  !  V1  =  j  ±  \  V 1  -  4r'~ 


(5.6) 


Thus,  waves  in  the  free  ocean  that  have  wavenumbers  in  the  range 


k  -  kRl  <  2kRVA0/c 


(5.7) 


(where  kR  is  the  wavenumber  of  the  resonance  [cg(kR)  =  c])will  be  reflected 
(i.e.,  reverse  their  group  velocity  as  seen  in  the  moving  coordinate  system) 
by  the  surface  currents.  In  the  ocean-fixed  coordinate  system  waves  having  the., 
wavenumbers  do  not  in  general  reverse  direction  but  appear  to  interact  strongly 

with  the  current  region. 
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The  varlution  of  the  wavenumber  with  position  in  the  current  pattern 
can  be  followed  with  the  aid  of  Figure  10.  Here  we  have  plotted  the  value 
of  the  reduced  wavenumber  k  as  a  function  of  the  local  value  of  the  reduced 
current  (6)  for  various  fixed  frequencies  (n  ■  constant).  In  this  plot  the 
unperturbed  ocean  is  the  vertical  axis  3  ■  0  and  the  point  representing  a 
given  wave  packet  moves  along  one  of  the  curves  n  *  constant  as  the  wave 
packet  traverses  the  current  region.  Reflected  waves  are  those  which  return 
to  the  6=0  axis  at  a  value  of  k  different  from  the  initial  value. 

For  example,  according  to  Figure  10,  waves  having  an  unperturbed  (reduced) 
frequency  n  “  0.24  can  be  reflected  if  3  dips  below  -  0.04.  Viewed  from 
the  ocean-fixed  coordinate  system,  these  waves  change  their  wavenumbers 
a«  a  result  of  a  "collision"  with  a  turning  point: 


<2 


and  lose  or  gain  velocity 


The  mean  square  amplitude  is  unchanged  (in  this  one-dimensional  problem) ,  but 


the  mean  square  slope  .can, change  within  the  limits: 
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FIGURE  10.  REDUCED  WAVENUMBER  IC(=C2I</g)  AS  A  FUNCTION 
OF  THE  CURRENT  PARAMETER  (3  (=4u/c)  FOR  FIXpD 
VALUES  OF  THE  FREQUENCY  ^(^CO/g) 
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These  changes  are  the  "residual"  changes  that  remain  after  the 
wave  has  ceased  Interacting  with  the  current  pattern.  The  changes  that 
take  place  during  the  course  of  the  Interaction  can  differ  considerably 
from  these  values.  However,  even  these  residual  values  can  be  significant. 
Consider  a  long  wavelength  Internal  wave  propagagatlng  at  0.5  meters/second 
along  a  thermocline  at  50  meters  depth.  An  amplitude  of  the  Internal 
wave  of  one  meter  would  create  a  surface  current  that  oscillates  between 
t  1.0  centimeters/second.  This  current  would  (In  the  absence  of  non-llnear 
effects)  he  capable  of  strongly  Interacting  with  surface  waves  having  wave¬ 
lengths  between  50  and  75  cm.  The  most  strongly-altered  waves  would  have 
their  wavelengths  changed  by  up  to  401.  their  group  velocity  by  up  to  20%, 
and  their  mean  square  slopes  by  up  to  80%. 


As  indicated  above,  the  moving  current  pattern  associated  with  a 
travelling  internal  wave  modifies  the  spectral  distribution  of  surface 
waves  in  two  ways.  Within  the  region  of  finite  current  gradient,  the  wave 
amplitudes  and  slopes  are  altered  according  to  Equations  (3.49)  and  (3.50). 

addition  to  this  local  effect,  some  of  the  waves  intercepted  by  the 
current  pattern  are  modified  by  the  scattering  process  described  in  Section  3 


and  leave  the  region  of  interaction  with  altered  wave  speeds.  As  these 
reflected  waves  propagate  away  from  the  current  pattern,  they  alter  the 


spectral  distributions  at  distances  remote  from  the  scattering  region.  To 
demonstrate  these  effects,  we  consider  the  effect  of  the  passage  of  a  simple 
gauss ian-shaped  internal  wave  packet  of  peak  amplitude  a  travelling  along 
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a  sharply-defined  thermocline  at  depth  D  (Figure  11).  ».  surface 

current  has  the  form 

MJ  _  u  a 

Cq  3  "  '  d  exp  <x  "  Cot)2/A2)  # 


and  we  evaluate  the  steady  state  spectral  distribution  at  differing 
spatial  locations  relative  to  the  current  packet. 

In  the  traae  of  the  Internal  wave,  R-w.ves  lmpi„ge  £ron  bo£h  directlons 
"<+>"  waves,  corresponding  to  wave  packets  moving  l„  the  positive  direction 

“d  h“Vl"8  8r°U|’  Vel0Cltl“  neater  than  that  of  the  Internal  wave,  and  »<-)•■ 
waves,  corresponding  to  those  .ovlng  slower  than  the  Internal  wave. 


In  Figure  12  we  show  x  - 
that  are  reflected  by  the 
these  plots. 


t  plots  of  wave  packet  trajectories.  The  waves 
current  regions  appear  as  discrete  bundles  in 


TO  demonstrate  the  modification  of  a  spectrum  of  linear  waves  by 
a  travelling  current  pattern,  we  assume  that  the  incident  waves  have  a 

spectrum.  This  spectral  for.  Is  chosen  to  simulate  that  of  a  saturated 
sea,  although  it  is  recognised  that  this  linear  wave  model  is  not  consistent 
with  the  concept  of  a  sea  limited  by  „o„-li„ear  effects.  In  Flg„re  13 
we  show  the  calculated  power  spectrum  (multiplied  by  kh  at  three 
different  points:  one  ahead  of.  one  behind,  and  one  directly  over  the 
internal  wave.  These  results  were  obtained  using  the  WKB  approximation 
formulas  (Equation.  4.21  and  4.22)  with  the  standing  wave  pattern  created 
by  the  reflected  wave  smeared  out.  The  data  ate  presented  In  terms 
of  dimensionless  variables  and  can  be  scaled  to  various  conditions.  For 
example,  for  an  Internal  wave  velocity  of  50  c./sec,  the  resonant 
wavenumber,  frequency,  and  wavelength  are:  k  -  0.097  cm'1,  „  .  9.8  se£-l  _ 
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Figure  12a.  Space-Time  Trajectories  for  the  Strongly 
Interacting  Gravity  Waves  Ol>(+)  waves) 
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Figure  12b. 


Space-Time  Trajectories  for  the  Stroup/  v 
Interacting  Gravity  Waves  (*(-)  waves) 
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Figure  12c.  Space-Time  Trajectories  for  the  Strongly 
Interacting  Gravity  Waves  (<f»^  waves) 
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IX> 


Figure  12d. 


Space-Time  Trajectories  for  the  Strongly 
Interacting  Gravity  Waves  (<*>(“)  waves) 
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Figure  13c.  Surface  Wave  Height  Spectrum  (times  k^)  Ahead 
of  the  Internal  Wave  (eMW  =  -  0.04) 


A  -  64  cm  (corresponding  to  the  dimensionless  values  ic  *=  0.25  and 

n  "  °*25-)*  For  the  sample  chosen  (0  -  -  0.04),  the  surface  current 
varies  from  zero  to  -  2.0  cm/sec  over  the  crest  of  the  wave. 

The  calculation  shows  a  number  of  interesting  effects.  First, 
over  the  internal  wave  crest  (Figure  13b),  the  spectrum  is  distorted  over  a  rather 
broad  range  of  wavenumbers.  The  broadening  is  less  ahead  or  behind  the 
internal  wave  (where  it  is  consistent  with  approximate  estimat  es  given  earlier). 

Another  effect  noticeable  is  the  lack  of  intense  caustic  formation. 

Although  a  single  wave  incident  on  the  pattern  would  show  a  large  amplitude 
enhancement  (infinite  in  this  approximation)  in  the  neighborhood  of  its 
turning  point,  this  is  not  true  when  the  incident  waves  are  spectrally 
distributed.  Here,  because  different  wave  packets  have  turning  points 

at  different  locations,  the  net  effect  on  the  spectrum  is  a  finite  dis¬ 
tortion  (see  Figure  13b). 

Outside  the  current  region  we  see  the  effect  of  reflected  waves. 

In  front  of  the  internal  wave  the  net  effect  is  an  apparent  depletion  of 
the  high  wavenumbers  (a  "smoothing").  This  results  from  the  fact  that, 
as  the  surface  waves  are  overtaken  by  the  internal  wave,  those  that  are 
"blocked"  by  the  current  are  accelerated  and  move  ahead  of  the  current 
pattern  with  decreased  wavenumber  (increased,  wavelength) . 

Because  the  spectrum  falls  cff  rapidly  with  lacreaslug  wavenumber,  the 
resulting  depletion  at  high  wavenumbers  is  much  more  noticeable  than  the 
corresponding  increase  at  low  wavenumbers. 

The  opposite  effect  occurs  astern  of  the  internal  wave.  Here  the 
reflection  process  increases  the  wavenumber  resulting  in  a  strong  increase 
in  the  high  wavenumber  part  of  the  spectrum  as  compared  to  a  fracti'onally- 
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weaker  decrease  in  this  low  wavenumber  parti  Thus  this  linear  model  would 
predict  an  apparent  high  wavenumber  smoothing  ahead  of  the  internal  wave 
crest  and  a  roughening  astema  Directly  over  the  crest  both  effects 
occur:  low  wavenumbers  are  enhanced,  high  wavenumbers  are  suppressed. 

The  morphology  discussed  above  applies  to  a  linear  sea.  Non— linear 
effects  in  the  real  ocean  can  be  expected  to  severely  limit  the  magnitude 
of  the  effects  induced;  however,  the  morphology  should  be  similar.  That 
is,  a  region  of  advancing  positive  current  gradient  (such  as  that  over 
the  bow  of  an  internal  wave  crest)  tends  to  accelerate  the  near  resonant 
waves  and  results  in  a  depletion  of  the  high  wavenumber  portion  of  the 
spectrum  (a  "smoothing").  Ar.  advancing  region  of  negative  current  gradient 
(such  as  in  the  convergence  zone  trailing  an  internal  wave  crest)  leads  to 
enhancement  of  high  wavenumbers  (a  "roughening"). 

One  can  infer  some  of  the  real  ocean  effects  by  following  the  motion 
of  an  individual  wave  packet  as  it  traverses  the  variable  current  region. 

In  Figure  14  we  show  histories  of  a  wave  packet  as  it  enters  the  current 
pattern  and  is  reflected.  The  mean  square  amplitude  enhancement  and  mean 
square  slope  enhancement  are  plotted  as  a  function  of  time.  In  a  truly 
linear  sea,  the  amplitude  and  slope  enhancement  depend  solely  on  the  local 
value  of  the  current  and  is  maximum  when  the  group  velocity  relative  to 
the  internal  wave  vanishes. 

The  enhancements  shown  in  Figure  14  are  much  larger  than  can  be 
expected  to  occur  in  a  fully-developed  sea,  however.  A  rough  estimate  of 


78 


DISTANCE  (x/A)  t  ENHANCEMENT  FACTOR 
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Figure  14.  Enhancement  Factors  and  Trajectory  for  a 

Wave  Packet  Having  an  Unperturbed  Wavelength 
of  95  cm  (the  resonant  wavelength  is  64  cm), 
Distance  is  measured  from  the  internal  wave 

CREST.  79 


the  reaJ  ocean  effects  to  be  expected  can  be  obtained  by  introducing  the 
concept  of  a  coherence  time.  This  is  the  time  required  for  wind,  wave-wave, 
or  dissipative  interactions  to  alter  the  wave  under  consideration.  Rough 
estimates  of  this  time  (measured  in  wave  periods)  are  shown  in  Figure  15 
and  indicate  that  coherence  times  in  the  range  of  10  to  100  wave  periods 
are  typical  for  a  saturated  sea.  The  wind  interaction  time  was  estimated 
from  Fhillips  (1966).  Two  estimates  are  given  for  the  non-linear  wave  in¬ 
teraction  time;  the  upper  one  according  to  Benjamin  (1967)  [t/T  =  k2a2/ir 

2  2 

and  k  a  ~  0.006  for  a  saturated  sea]  and  the  lower  one  according  to  Has- 

2  2 

selman  (1967)  [t/T  ~  with  the  mean  square  slope  fitted  to  the  data 
of  Cox  and  Munk] . 

For  example,  consider  the  gausslan-shaped  internal  wave  discussed 
previously.  This  wave  had  a  half  width  A  *  50  meters,  a  peak  amplitude 
a  =  2  meters,  and  was  assumed  to  be  travelling  at  50  cm/sec  along  a  sharp 
thermocline  at  50  meters  depth.  We  follow  a  surface  wave  which,  in  the 
unperturbed  ocean  has  a  wave  length  equal  to  95  centimeters  (wave  period 
of  0.78  sec.)  and  which  travels  at  a  group  velocity  equal  to  60.8  cm/sec 
(i.e. ,  it  is  overtaking  the  internal  wave  at  an  initial  relative  velocity 
of  10.8  cm/sec). 

According  to  Figure  14  this  wave  will  be  decelerated  to  the  internal 
wave  speed  by  the  time  it  has  penetrated  to  a  point  24  meters  behind  the 
internal  wave  crest.  For  such  a  wave  the  width  of  the  "full  wave"  region 
(as  defined  by  the  point  where  the  Airy  integral  argument  c,  is  equal  to 
unity  -  see  Eq.  3.59)  is  given  by  the  expression 
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U*  {cm./sec.) 

Figure  15.  estimated  interaction  times  for  various  effects, 

HERE  U*  IS  THE  TURBULENT  FRICTION  VELOCITY  AND 
t  THE  WAVE  PERIOD. 
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Ax  *  j4c/(k^du/dx)  j 

where  du/dx  is  the  local  current  gradient.  At  the  reflection  point  in 
our  example,  the  current  gradient  is  about  -  4  x  10  4  sec  and  Ax  is 

about  75  cm. 


If  the  sea  were  truly  described  by  linear  theory,  the  mean/square 
amplitude  and  slope  would  have  been  enhanced  by  the  factors  6  and  10, 
respectively,  at  the  edge  of  the  "full  wave"  region.  However,  the  time 
for  the  wave  packet  to  travel  from  a  point  at  the  edge  of  the  current 
region  (say  x  *  -  50  meters)  to  the  point  of  maximum  enhancement  (the 
edge  of  the  "full  wave"  region)  is  very  long:  At  sss  6  A/c  =  600  seconds 
or  about  770  wave  periods.  According  to  Figure  15,  it  is  unlikely  that 
such  a  wave  will  have  such  a  long  coherence  time  except  when  the  wind 
speed  is  quite  low.  Thus  it  is  more  appropriate  to  ask  what  happens  to 
a  wave  which  is  born  within  the  current  pattern  and  which  travels  a 
relatively  short  distance  from  its  birth  place. 

Consider  a  wave  generated  at  the  point  A  in  Figure  14  (26  meters 
beind  the  crest  of  the  internal  wave).  Such  a  wave  has  a  wavelength  of 
74.2  cm  and  a  group  velocity  of  52.2  cm/sec  (2.25  cm/se_  relative  to  the 
internal  wave).  In  50  seconds  (64  wave  periods),  it  will  move  co  a  point 
1.12  meters  closer  to  the  internal  wave  crest.  During  this  time  its  mean 
square  amplitude  and  slope  will  be  enhanced  by  the  factors  1.23  and  1.30 
respectively.  These  values  are  probably  more  realistic  for  the  normal 
state  of  the  ocean  than  the  values  quoted  in  the  previous  paragraph. 
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It  is  to  be  noted  that  the  enhancement  effects  discussed  above  are 
only  one  of  two  effects  that  result  in  modification  of  ocean  wave  spectra. 
Even  in  the  absence  of  enhanced  amplitudes,  a  spectrum  will  be  altered 
because  of  che  simultaneous  wavenumber  changes  induced  by  the  current  gra¬ 
dient  (refraction).  In  order  to  assess  the  simultaneous  effects  of  these 
two  phenomena,  as  well  as  to  account  for  the  finite  coherence  time  effects 
and  the  effects  of  non-linear  saturation,  it  is  necessary  to  develop  a 
formalism  for  treating  the  spectral  deformation  in  terms  of  a  spectral  re¬ 
laxation  or  transport  equation.  The  development  of  such  a  model  will  be 
the  subject  of  subsequent  reports. 

The  spectral  modifications  predicted  by  the  present  model  are  summarized 
Figure  16.  The  phasing  of  the  "smoothed"  region  (i.e.,  the  mechanical 
slicks")  in  this  linear  wave  model  is  just  opposite  to  that  one  would  ex¬ 
pect  as  a  result  of  enhanced  viscous  damping  in  the  convergence  regions  due 
to  an  accumulation  of  organic  material  or  other  debris. 

Although  one  is  tempted  to  draw  conclusions  concerning  the  phase  re¬ 
lationship  of  regions  of  roughening  and  smoothing  in  the  ocean  based  on 
this  linear  model,  it  must  be  remembered  that  the  actual  predicted  effect 
is  one  of  spectral  distortion  with  simultaneous  enhancement  of  some  wave¬ 
lengths  and  suppression  of  others.  The  non-linear  effects  present  in  real 
saturated  ocean  spectra  are  not  likely  to  affect  all  wavenumbers  equally. 

Thus  the  model  needs  to  be  extended  to  include  the  effects  of  wind  genera¬ 
tion  of  waves  and  of  non-linear  wave-wave  interactions,  including,  wave 
breaking,  before  one  can  make  confident  predictions  of  real  ocean  wave 
spectra  modifications. 
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Figure  16,  Morphology  of  Mechanical  Wave  Modification 
in  a  Linear  Sea 


APPENDIX:  MICRO-STRUCTURE  OF  THE  TURNING  POINT  AND  APPLICABILITY  OF  THE 


WKB  APPROXIMATION 


As  pointed  out  In  the  WKB  approximation  of  Section  4,  there  exists 
the  possibility  of  reflected  waves  at  the  WKB  turning  points,  i.e, ,  the 
points  where  the  wave  amplitude  becomes  infinite.  According  to  the  dispersion 
relation  [Equations  (3.44)  and  (3.45)],  the  reduced  wavenumber  k  has  an 
imaginary  part  (corresponding  to  evanescent  waves)  when 


1  -  4n(l  -  0)  +  4nV  <  0 


(A.  1) 


in  other  words,  when  the  reduced  frequency  n  lies  between  the  critical  values: 


and 


1  -g 
202 


"2  '  ^ 

202 


1  -  /l-02/(1-  0)2 


I - - - 

1  +  jl-  0V(1-  0)2 


(A.  2) 


(A.  3) 


that  is,  when 


"l  S  ”  S 


(A. 4) 


the  wave  cannot  propagate.  Thus,  for  small  values  of  g,  evanescent  waves  are  ex¬ 
pected  when 


1  *  n  1 
“7  d  n  >  — — 

02  4 


(A. 5) 


which  is  obtained  by  expanding  Equation  (A, 2).  Of  primary  interest  to  us  is 
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the  low  frequency  limit,  since  we  learned  earlier  that  at  this  limit  the 
reflected  wave  amplitude  becomes  infinite  in  the  WKB  approximation. 

To  expand  on  this  result,  let  us  consider  the  group  velocity  of 
the  surface  wave  in  this  (moving)  coordinate  system.  Using  the  dispersion 

relation  (Equation  (5.2)]  for  the  reduced  frtequency,  and  neglecting  terms 
gradratic  in  3> 


n  *  -  < 


(l  -0)  ±  \Zk~ 


(A.6) 


we  obtain  for  the  group  velocity  (relative  to  the  internal  wave  velocity  c) 

^  =  dn  =  _  _  j_ 

c  '  2/r  (A-7‘> 

where  the  upper  sign  corresponds  to  the  low  frequency  mode  of  interest;  thus 


=  -  1  +  ~~  +  e  (a.  7b) 

/4k 

and,  at  the  point  n  -  n1~l/[4(l-0)]  where  =  1/  [4(1-20)  ],  the  group  velocit 
is 

c  n 

-  -  ( 1  -  0 )  +  Vl-  2  0  a  0  .  (A.  8) 

The  point  n  -  is  thus  a  turning  point  and  corresponds  to  specular  or  total  re¬ 
flection  of  the  gravity  wave  by  the  surface  current  at  points  where  g  =  i  _  1/4 
[see  Equation  (4.23)]. 


In  the  evanescent  region  the  wavenumber  is  complex  and  may  be  written  In  the 

form 


K 


(A.  9) 
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which  when  substituted  Into  the  dispersion  rsl.tlon  Equation  <3. 11)  yields 

values  for  tc  and  k  : 

K  I 


K  -  1-2(1-  g) 

R  2(1-  28)  (A. 10) 

and 

2  (1-  g)n  -  i 

1  (1  -  20) 2  ’  (A. 11) 


For  small  8  these  expressions  sre  [  In  term,  of  the  critical  frequency 
1,1  "  4(1-6)  |  ! 


and 


(A. 12) 


(A. 13) 


A  measure  of  the  applicability  of  the  WKB  approximation  is  the  relative  magni¬ 
tude  of  the  "skin  depth"  or  penetration  depth  at  the  turning  point  compared  to 
the  dimension  of  the  current  pattern.  The  penetration  length  is  of  order 


It  is  convenient  to  evaluate  this  length  near  the  resonant  condition,  i.e., 
near  the  value  of  kr  for  which  I  /£7k*  -  c.  For  an  internal  wave  propagating 
along  an  interface  at  depth  D  between  two  uniform  layers  that  differ  in  density 
by  an  amount  Ap,  we  may  write  (since  c2  -  g(Ap/p)D) 


k*  -  p/DAp 


(A. 15) 
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(A. 16) 


As  6  Increases,  higher  and  higher  frequencies  can  propagate.  Thus, 
the  negative  parts  of  the  current  will  reflect  the  lowest  frequencies  and 
the  more  negative  is  the  value  of  the  8  in  the  wave  trough,  the  stronger  will 
the  extinction  of  the  surface  waves  be.  The  value  that  (-8)  must  exceed  the 
critical  value  (l-4ri)  in  order  that  the  penetration  depth  be,  say,  5%  or  less 
of  the  thermocline  depth  D  is  given  by: 


-05%  =  1  -  4n  +  0.0016  (A. 17) 

where  we  have  taken  an  air-water  interface  (Ap/p  ■  10  ^).  Thus,  small  changes 
in  8  will  cause  a  wave  that  could  just  penetrate  the  entire  current  region  to 
be  reflected  in  quite  short  distances.  This  change  is  small  compared  to  typical 
variations  of  8  expected.*  Thus,  for  most  of  the  waves  we  can  assume  that  if 
they  are  reflected  at  all,  they  will  be  reflected  from  highly  localized  regions. 


*  The  internal  wave  height  Yj  is  given  by 

III!!  _  8(E) 

D  KD 

For  the  typical  case  KD  $  1,  significant  wave  heights  (Yj>  .01D)  are  always 
associated  with  8  variations  large  compared  to  that  in  Eq.  (A. 17). 
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